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Remote sensing has become important in the ocean sciences, especially for 
research involving large spatial scales. To estimate the in-water constituents through 
remote sensing, whether carried out by satellite or airplane, the signal emitted from 
beneath the sea surface, the so called water-leaving radiance (LJ), is of prime 
importance. The magnitude of L w depends on two terms: one is the intensity of the 
solar input, and the other is the reflectance of the in-water constituents. The ratio of 
the water-leaving radiance to the downwelling irradiance (E d ) above the sea surface 
(remote- sen sing reflectance, R rs ) is independent of the intensity of the irradiance input, 
and is largely a function of the optical properties of the in-water constituents. 

In this work, a model is developed to interpret R rs for ocean water in the 
visible-infrared range. In addition to terms for the radiance scattered from molecules 
and particles, the model includes terms that describe contributions from bottom 
reflectance, fluorescence of gelbstoff or colored dissolved organic matter (CDOM), 
and water Raman scattering. By using this model, the measured R rs of waters from the 
West Florida Shelf to the Mississippi River plume, which covered a [chi a] 
(concentration of chlorophyll a) range of 0.07 - 50 mg/m 3 , were well interpreted. The 
average percentage difference ( a.p.d .) between the measured and modeled R rs is 
3.4%, and, for the shallow waters, the model-required water depth is within 10% of 
the chart depth. 

Simple mathematical simulations for the phytoplankton pigment absorption 
coefficient (a*) are suggested for using with the R r5 model. The inverse problem of 
R rs , which is to analytically derive the in-water constituents from R rs data alone, 
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can be solved using the functions without prior knowledge of the in-water optical 
properties. More importantly, this method avoids problems associated with a need for 
knowledge of the shape and value of the chlorophyll-specific absorption coefficient. 
The simulation was tested for a wide range of water types, including waters from 
Monterey Bay, the West Florida Shelf and the Mississippi River plume. Using the 
simulation, the ^-derived in-water absorption coefficients were consistent with the 
values from in-water measurements (r 2 > 0.94, slope — 1.0). 

In the remote-sensing applications, a new approach is suggested for the 
estimation of primary production based on remote sensing. Using this approach, the 
calculated primary production (PP) values based upon remotely sensed data were very 
close to the measured values for the euphotic zone (r 2 = 0.95, slope 1.26, and 32% 
average difference), while traditional, pigment-based PP model provided values only 
one-third the size of the measured data. This indicates a potential to significantly 
improve the accuracy of the estimation of primary production based upon remote 
sensing. 
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INTRODUCTION 


Introduction to Remote Sensing 

Remote sensing has become an important tool in the study of ocean science, 
especially for the synoptic research of large bodies of water, because a large region of 
the world oceans can be viewed by a remote sensor in a very short time. For 
example, the Nimbus-7 and TIROS satellites had repeat cycles of 2 days at the 
equator with > 1 km pixel diameter; at the other end of the scale, however, a blimp 
can transit Tampa Bay in mere hours measuring at high resolution ( — 8 x 8 m 2 ) and 
lingering long enough to observe changes with time. 

The main purpose of oceanic remote sensing in the visible-infrared region of the 
spectrum is to remotely derive concentrations of certain constituents in the water 
(principally phytoplankton, dissolved organic matter and suspended inorganic 
sediments), to estimate the oceanic primary production and to understand ocean 
circulation. The derivation of these constituents through remote sensing, whether by 
satellite or aircraft, depends on the signal emitted from beneath the sea surface, the so 
called water-leaving radiance L w (\) (definitions of variables are in the List of 
Symbols; the functionality of one variable to another may not be explicitly expressed 
unless it is needed for further clarity), or the spectrally digitized water color. There 
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are two approaches to formulate this derivation from measured L „ : one is empirical, 
and the other is analytical. The empirical approach, such as the Coastal Zone Color 
Scanner (CZCS) algorithm for pigments [Gordon et al. 1983], is most often used at 
present. 

The empirical approach is based on observations of the change of water color 
(spectral ratios) as change of water constituents, and statistical methods are used to 
find the apparent relationship of one variable to another. Ocean waters show a wide 
variation of colors, ranging from dark blue through blue/green, green, and 
green/yellow [Jerlov 1976], For clear sky days, blue water is usually found in the 
open ocean, where the water is clear, deep, and very few light absorbers such as 
phytoplankton are present. Pure sea water (which means the sea water without 
suspended particles and dissolved organic materials) is a good absorber for 
wavelengths longer than 570 nm [Smith and Baker 1981]. For upwelling regions, 
water color turns to green due to the increased concentration of phytoplankton and its 
associated degradation products (gelbstoff and detritus), which are strong absorbers in 
the blue/green bands [Jeffrey 1980, Bricaud et al. 1981, Roesler et al. 1989]. For 
coastal (especially shelf) regions, gelbstoff (yellow substance) or colored dissolved 
organic matter (CDOM) is abundant due to land run-off and absorbs strongly in the 
blue, turning the water to greenish yellow colors. 

Based on the above observations, it is easy to qualitatively scale the 
concentrations of in-water constituents such as chlorophyll a ([chi a]) and gelbstoff 
([YS]) based on ocean color, and this is the basis for the CZCS algorithm. 
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The CZCS algorithm is a power law of the ratio of L (V (443)/L W (550), which was 
empirically developed. Its accuracy in deriving pigment concentration is about a factor 
of 2 [Gordon and Morel 1983] for the oligotrophic ocean, and is much worse for 
coastal waters. This is because the water color is affected by constituents in addition 
to chlorophyll alone, especially in coastal areas. Water color signals, as observed 
immediately above the ocean’s surface, depend upon the incident irradiance and the 
concentration and distribution of the in-water constituents. The waters will appear 
bluer on a bright, sunny day than on an over-cast cloudy day as there is relatively 
more blue light entering the sensor from reflected sky light rather than cloud light. At 
the same time, in-water constituents which can affect the water color include 
chlorophyll, gelbstoff, detritus, bacteria, and inorganic particles. Transpectral 
responses (water Raman scattering, gelbstoff fluorescence) and bottom reflectance 
further complicate matters. This is why simple ratios of L w ’s do not always work well 
to derive [chi a] [Carder et al. 1986], especially for coastal waters, when there may 
be a strong bottom influence. Figure 1 illustrates the schematic relationships between 
water-leaving radiance and in-water components. 

Since there are many components in the water which can affect the water color, 
it is necessary to separate their individual influences if we want to improve the 
algorithm. This favors the analytical approach. For the purpose of deriving in-water 
constituents, first it is essential to understand how each component affects the optical 
signal. Only after this step, can we have the potential to analytically derive the 
individual component from the remotely measured signal alone. So, there are two 
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chi a , YS, ... 



Figure 1 . Schematic Relationships between Water Color and In-water Constituents. 

steps to the analytical approach as pointed out by Morel [1980]. For the remote- 
sensing technique, first is the forward step (or the forward problem), i.e., to interpret 
the observed color based on the measurement of each component (solid arrow in the 
box of Figure 1); second is the backward step (or the inverse problem), i.e., to derive 
each component just from the remotely measured signal (dashed arrow in the box of 
Figure 1). 

The forward problem can largely be solved based upon measurements of the 
optical components. If a few parameters prove difficult to measure, they may be 
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simulated based on theoretical grounds. This step is important not only for the 
understanding of the optical scenario, but also for the guidance of future theoretical 
and experimental studies. 

In the forward problem, as upwelling light depends upon the intensity of input 
irradiance, subsurface irradiance reflectance (K(O-)) and remote-sensing reflectance 
(/?„) are defined as the ratio of upwelling irradiance to downwelling irradiance just 
below the surface and the ratio of the water-leaving radiance (L w .) to the downwelling 
irradiance (E d ) above the sea surface, respectively. In this way, R( O') and R rs will be 
independent of the intensity fluctuation of the input irradiance, and largely a function 
of the optical properties of the in-water constituents. Since E d can be very precisely 
predicted from models [Gregg and Carder 1990, Bishop and Rossow 1991], it is easy 

to calculate L w if we know how to accurately determine R„. 

Optical models have been developed for the subsurface irradiance reflectance 
*( 0) [Gordon et al. 1975, Morel and Prieur 1977, Gordon 1989b, Kirk 1984, 1991], 
but satellites measure the upwelling radiance above the water surface. To use satellite 
radiance data, it is necessary to know how to interpret remote-sensing reflectance 
instead of irradiance reflectance. Models were suggested by Carder and Steward 
[1985], and Peacock et al. [1990] to explain the measured R rs , but in each case no 
bottom reflectance, gelbstoff fluorescence, or water Raman scattering were included. 
Also, in these works a somewhat arbitrary "Q factor" [Austin 1974] was used. 

For the open ocean, in-water constituents are dominated by phytoplankton 
particles and their derivatives. As a result, R rs is easier to interpret, and is largely 
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dependent on the concentration of chlorophyll. For these "case 1" [Morel and Prieur 
1977] waters, Monte Carlo simulations [Gordon et al. 1988, Morel and Gentili 1993] 
have been made for R n (Q ). However, it is not clear whether these simulations can be 
used for "case 2" waters due to the different volume scattering functions and larger 
chlorophyll a concentrations. 

For coastal waters, as discussed previously, components such as land run-off, 
suspended sediments and bottom reflectance, which do not co-vary with [chi a], 
confound the data interpretation. Also, the concentration of gelbstoff in coastal waters 
is frequently high and can fluoresce strongly in the blue-green region [Coble et al. 
1990 and references cited], but it is not clear how this fluorescence might affect the 
remote-sensing reflectance. Further, it has been found that water Raman scattering 
can contribute as much as about 10% to the sub-surface irradiance reflectance in the 
open ocean [Stavn 1988, 1990], but how this would influence the R rs is not clear. 

For shallow coastal waters, part of the water-leaving radiance is frequently due 
to bottom reflectance, so, for the derivation of [chi a ] from R rs , the influence of the 
bottom reflectance needs to be removed first. How we express the bottom influence in 
R„ needs further development. 

For water-depth measurements [Polcyn et al. 1970, Clark et al. 1987] or 
bottom-feature mapping [Lyzenga 1978], the diffuse attenuation coefficient [Clark et 
al. 1987] or an "effective" attenuation coefficient [Lyzenga 1978] was used. For 
shallow water R( 0‘), Spitzer and Dirks [1987] used 3.7 a (the total absorption 
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coefficient) to represent the "effective" attenuation coefficient. It is not clear, 
however, how to express the "effective" attenuation coefficient for R rs . 

The inverse problem is to use remotely measured data to derive the absorption 
and scattering properties and to estimate in-water constituents such as [chi a], [YS], 
etc. These derived values can be used to model the spectral light at depth 
[Sathyendranath and Platt 1988] and/or to estimate primary production [Morel 1988, 
Balch et al. 1992]. Analytically, not much improvement has been made for the 
inverse problem in the past studies, mainly because it is not easy to accurately express 
the change of each component by a few parameters. 

In the application of remote sensing, empirical and semi-analytical approaches 
have been studied in the estimation of [chi a], and efforts have been made to estimate 
the primary production based on pigment concentration and empirical chlorophyll- 
specific absorption coefficients [Morel 1991, Platt et al. 1991, Balch et al. 1992]. 

But, as the chlorophyll-specific absorption coefficient varies significantly from sample 
to sample [Morel and Bricaud 1981, Bricaud and Stramski 1990, Carder et al. 1991, 
Carder et al. 1994], the concentration of chlorophyll cannot be derived accurately 
from remote sensing without prior knowledge of the chlorophyll-specific absorption 
coefficient, which will further influence the primary production. Therefore, for cases 
where only water-leaving radiance (or remote-sensing reflectance) is available, better 
methods are needed for the derivation of the in-water absorption coefficient and the 
estimation of primary production. 
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The Present Study and Structure 

This contribution uses in-water and remote-sensing optical data to develop 
methodologies that can be used with present and new aircraft and spacecraft sensors 
to quantify biological, chemical and physical properties of the surface waters from 
open ocean to coastal waters. 

In this study, a mathematical model is first developed to interpret the remote- 
sensing reflectance (solid arrow in the box of Figure 1) based on previous theoretical 
and experimental studies. In the model, terms to describe contributions from bottom 
reflectance, gelbstoff fluorescence and water Raman scattering are developed. In 
addition, hyperspectral R n measurements taken from waters ranging from the West 
Florida Shelf to the Mississippi River plume are selected to demonstrate the model. 
Differences between the modeled and measured R rs are discussed. 

Second, new approaches are suggested for the inverse problem (dashed arrow in 
the box of Figure 1), which include mathematical and empirical simulations for the 
absorption curves of phytoplankton pigments (a^). With these simulations, R„ or 
upwelling radiance above the surface (Z,„(0 + )) can be modeled when there are no in- 
water measurements available. Alternatively the in-water absorption coefficient can be 
derived from remote signals alone, and the spectral light field at depth as well as the 
primary production can be calculated. 

Part A of this manuscript describes the optical basis of remote sensing and the 
model development. It is divided into 4 chapters: Chapter 1 , the optical properties of 



9 


Table 1. Cruise Information where Data Were Collected and Used in this Study. 


Cruise Name 

Location 

Period 

Water Type 

Monterey (MA) 

Monterey Bay 

Sep. 1989 

green 

Monterey (MB) 

Monterey Bay 

Oct. 1989 

green 

TBX1 (TA) 

West Florida Shelf 


shallow, green 

ML-ML1 (ML) 

North Atlantic 

May-June 1991 

blue 

TBX2 (TB) 

West Florida Shelf 

May 1992 

green, blue 

GOMEX (GO) 

Gulf of Mexico 

April 1993 

green, blue 

COLOR (CO) 

Gulf of Mexico 

May-June 1993 

dark, green 

SHAMU (SH) 

Tampa Bay 

Jan. 1993 

shallow 


Note: Characters in parentheses will be used in text and tables as a short name. 


pure sea water, gelbstoff and phytoplankton pigments, and the relationships between 
the apparent and inherent optical properties; Chapter 2, model development for 
remote-sensing reflectance; Chapter 3, data used and measurement methods; Chapter 
4, model results; and Chapter 5, conclusions of part A. Part B presents the ways to 
solve the inverse problem and discusses some applications of remote sensing: 

Chapter 6, the methods suggested to solve the inverse problem; Chapter 7, 
comparisons of the derived absorption coefficient from R rs inversion and that from K d , 
Chapter 8, discussion of the remote-sensing applications, with a focus on primary 
production; Chapter 9, conclusions of part B; and Chapter 10, summary of this study 
and expectations for future work. 

Data used in this study came from cruises during the period of 1989 to 1993, 
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which covered a wide range of water types, such as high-latitude North Atlantic 
waters; sub-tropic waters of Monterey Bay; gelbstoff rich, shallow waters of the West 
Florida Shelf; clear waters in the Loop Current; and the phytoplankton bloom waters 
in the Mississippi River plume. Table 1 summarizes the cruise information, Figure 2a 
shows the locations of the data collection, and Figure 2b zooms in on the stations 
where R n curves are used for the model demonstration. 



Figure 2b. Stations where Data Were Used to Demonstrate R n Model. 




PART A 


THE OPTICAL BASIS AND THE FORWARD PROBLEM 
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CHAPTER 1 

APPARENT AND INHERENT OPTICAL PROPERTIES 


1.1 Introduction 


There are two kinds of optical properties which relate to the in-water 
constituents [Presendorfer 1976]: apparent optical properties (AOPs) and inherent 
optical properties (IOPs). AOPs are functions of both the light field and the in-water 
constituents, whereas IOPs have no relation to the light field and are solely dependent 
on the distribution and optical character of the in-water constituents. AOPs are easy to 
measure routinely in optical ocean studies, but difficult to relate to in-water 
constituents. For IOPs, the situation is reversed. 

One goal of optical oceanography is to derive the IOPs from the AOPs, then to 
evaluate the physical, biological and/or chemical conditions of the water and global 
ocean. This process involves some optical quantities that can be measured directly, 
such as the radiance and irradiance fields. The connections between these optical 
quantities and the IOPs are through the AOPs. In this chapter, we will briefly discuss 
these optical quantities and properties. 
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1.2 Optical Quantities Describing the Light Field 

The optical quantities used most widely in describing the light field are radiant 
intensity, radiance, vector irradiance, and scalar irradiance: 

Radiant intensity, I(z,Q,<p), is defined as the flux of radiant energy per solid 
angle in a specified direction, and has the units W/sr. 

Radiance, L(z,Q,<p), is the radiant intensity per unit area, i.e., L = dl/dA. The 
area dA is perpendicular to the light beam. Radiance has the units W/sr/m 2 . 

Vector irradiance describes the total radiant intensity incident on a plane 
surface. In ocean optics, two useful vector irradiances are defined: downwelling 
irradiance, E d , and upwelling irradiance, E u . These two are the values of the 
irradiance on the upper and lower faces, respectively, of a horizontal surface. 

The relationships between radiance and the downwelling and upwelling 
irradiance are as follow: 


( 1 ) 


E d = J'^L(d,(p)cos(6)da) , 

E u = -f ^L(0,<p)cos(d)da> , 

where the first expression is an integral over the up-looking hemisphere. The negative 
sign before the second integral is used to force E u to be positive. Figures 3a and 3b 


show examples of E d and E u spectra versus depth, respectively. 

Scalar irradiance, E 0 , is the total radiant intensity coming from all directions, 




400 500 600 700 

wavelength (nm) 

(b) 

Figure 3. Examples of Downwelling (a) and Upwelling (b) Irradiance Spectra. 
(Courtesy of Dr. J. L. Mueller) 
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measured at a point in space, so 

E a = / L(d,(p)do) . (2) 

This term is appropriate for evaluating the flux received by a phytoplankton cell. 

It is useful to divide the scalar irradiance into downwelling (E^) and upwelling 
( E <J components, which are integrals of the radiance distribution over the upper or 
lower hemisphere, respectively, 

E od = L L(Q,<p)du , 

J r 2 * (3) 

E ou = J 2 m<p)dG> . 

Both vector and scalar irradiance have units of W/m 2 . 

1.3 Apparent Optical Properties (AOPs) 

AOPs come from simple manipulation of the optical quantities and are jointly 
dependent on the IOPs and on the distribution of the ambient light field. The 
important AOPs are the irradiance reflectance ( R ), remote-sensing reflectance (R rs ), 
and the diffuse attenuation coefficients for downwelling (KJ and upwelling (K u ) 
irradiance. 


1.3.1 Irradiance Reflectance, R 


Irradiance reflectance (or irradiance ratio), R(z), is the ratio of the upwelling to 




Figure 4. Examples of Irradiance Reflectance for Different Waters. 
(Courtesy of Dr. J. L. Mueller) 


the downwelling irradiance at a given depth z in the light field: 


R(z) = 


*.o 

EJ® ‘ 


R has no units. Figure 4 shows examples of sub-surface R for two different 


water types, where [chi a] was 0.12 mg/m 3 for GO10, and 7.4 mg/m 3 for GO04. 


1.3.2 Remote-sensing Reflectance, R n 


Remote-sensing reflectance, R rs , is generally defined as the ratio of the 
upwelling radiance to the downwelling irradiance at a given depth 
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*„(z,6 w ,<p) = 


LjLzfi^y) 

E/z) 


(5) 


where 9 W is the in-water zenith angle and <p is the azimuth angle from the solar 
plane, respectively, for the upwelling radiance L u . Figure 5 shows the light geometry. 
R n has the units of steradian' 1 (sr 1 ). Since a satellite or aircraft sensor measures the 
radiance leaving the water surface, a very useful form of R„ is 


p> 


^ w (Q%e a ,<p) 

Ejsn 


( 6 ) 


i.e., the ratio of the water-leaving radiance to the downwelling irradiance just above 
the surface. Where 0 a is the in-air zenith angle of the water-leaving radiance, and 
sin(0J = n w * sin(0 w ) (n„ is the refractive index of sea water). In most cases, 0 a is 
within 30° and <p is about 90° from the solar plane for ship-board measurements. In 
the following text, 0 and <p may not be explicitly shown for the sake of brevity. 
Figure 6 shows examples of R rs of different waters, with [chi a] of 16.2, 2.4, 0.6 and 
0.5 mg/m 3 for stations GOOl, GO05, GO20, and G022, respectively. 


1.3.3 Diffuse-attenuation Coefficient of Irradiance 


This property is divided into two components; one for the downwelling 
irradiance, K d , and another for the upwelling irradiance, K u . These represent the rate 
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of change of the natural log of the irradiance with depth, 

„ d UEJiz)] 

dz ( 7 ) 

d In [E u (z)] Kn 

K - ~~ — • 

The minus before the derivative is to force K d and K u to be positive, as z is positive 
downward. The units for K are nr 1 . 

Similar to Lambert-Beer’s Law [Gordon 1989a], this definition clearly indicates 
that, in any absorbing and scattering medium (in this case sea water), all the 
quantities describing the light field change with depth, z. The change is typically a 
decrease in E d and E u with increasing depth, in an approximately exponential manner. 
So it is convenient to specify the rate of change as the change of the natural logarithm 
of the value with depth. When the rate of change and value at one depth are known 
(e.g. K d and E d (Qr)), it is possible to estimate the value at another nearby depth for a 
well mixed ocean. 


1.3.4 Distribution Function and Average Cosine 

Distribution functions (D d (z) and D a (z)) and average cosines (jl^z) and ~ix u (z)) 
for downwelling and upwelling light, respectively, are defined as 
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D/z) = - 


M z > 


EJz) 

E/z)' 


D u (z) « 


1 

V u (z) 


E ou (z) 
E u (z) ' 


( 8 ) 


These factors are simple means of characterizing the shape of the radiance 
distribution. For example, for completely diffused upwelling light, D u = 2.0, n u = 
0.5. For a collimated downwelling light beam, D d = n d = 1. 

In most oceanic cases, both D d and D u increase with depth, and DJD d « 2 
[Gordon et al. 1975]. For the downwelling distribution function just beneath the 
surface, D/0) * l/cos(/j, with an error of less than 3% for a clear-sky situation 
[Gordon 1989a], where j is the sub-surface solar zenith angle. 


1.4 Inherent Optical Properties (IOPs) 

IOPs are those that are independent of the ambient light field (i.e. independent 
of the intensity and the angular distribution). The important IOPs of any medium are 
the absorption coefficient, a, the scattering coefficient, b, and the volume scattering 
function (VSF), j8(a), which expresses the magnitude and angular distribution of 
scattered photons. The angle a describes the direction change of the scattered photons 
from the input photons. An additional inherent optical property is the beam 
attenuation coefficient, c, given as 

c = a + b. 

The beam attenuation coefficient thus represents the fractional loss of photons from a 
collimated light beam due to absorption and scattering, with units m' 1 . 
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wavelength (nm) 


Figure 7a. Absorption Coefficients of Pure Water. 


The absorption and scattering coefficients of water as a whole, at a given 
wavelength, are equal to the sum of the individual coefficients of the components 
present [Gordon et al. 1980]. 


1.4.1 Absorption Coefficient 


There are three major components for the absorption coefficient of ocean water: 
water itself, a w \ suspended particles, a p \ and gelbstoff, a g . a p can be further divided 
into two components: the absorption coefficient of phytoplankton pigments, a ^ and 
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the absorption coefficient of detritus, a d . So: 

a = a w + a t + a p or a = a w + a g + a d + a*. 

14.1.1 Absorption Coefficient of Water Itself . The absorption coefficient of water 
itself is one of the most important optical components. Figure 7a shows three 
measured absorption spectra. Water absorbs light very weakly in the blue and green 
regions of the spectrum, but increases markedly as the wavelength rises above 570 
nm, becoming quite significant in the red region. Also notice in the figure that there 
are some differences in a w among researchers. The main reason is the difficulty in 
obtaining "pure" water for experiments (also see the discussion in Smith and Baker 
(1981)) and the long path length needed for accuracy at blue wavelengths. We will 
use the a w values from Smith and Baker [1981] as their work is more recent and 
widely used. But, it must be remembered that the accuracy of their a K was only 

+25% to -15% between 380 to 800 nm. 

For wavelengths in the visible, a w does not change much between fresh and salt 
waters [Morel 1974]. However, due to the dissolved ions in the sea water, a w can 
change markedly for wavelengths in the UV range [Morel 1974], 

Hojerslev and Trabjerg [1990] reported that a w varies with temperature. 
However, their results differed from Pegau and Zaneveld [1993]. So, questions about 
how a w (k) varies with temperature remain. 
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wavelength (nm) 

Figure 7b. Examples of Absorption Coefficients of Gelbstoff. 

1-4. 1.2 Absorption Coefficient of Gelbstoff . Figure 7b shows examples of the spectral 
absorption coefficient of gelbstoff. It is a strong blue absorber with exponential 
decrease with wavelength. This absorption coefficient can be expressed as [Bricaud et 
al. 1981, Carder et al. 1989, Roesler et al. 1989] 

a g (X) = a g (440)e '■V x ' 440) . (9) 

The spectral slope S g varies with different dissolved organic materials, which also 
varies slightly with wavelength range [Krijgsman 1994]. At 440 nm, Carder et al. 
[1989] found that S g « 0.011 nm' 1 for a marine humic acid and 0.019 nm' 1 for a 
marine fulvic acid. The average of S g for ocean waters is reported as 0.014 nm' 1 
[Bricaud et al. 1981, Kishino et al. 1984], 
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wavelength (nm) 

Figure 7c. Examples of Absorption Coefficients of Phytoplankton Pigments. 

1.4. 1.3 Absorption Coefficient of Phytoplankton Pigments . Figure 7c presents 
examples of the spectral absorption coefficient of phytoplankton pigments, which are 
the basis of oceanic photosynthesis. This spectrum typically has two peaks, one 
around 440 nm, another around 675 nm. These are due to the presence of chlorophyll 
a. The shoulders above 440 nm are due to the presence of the accessory pigments, 
such as chlorophyll b , chlorophyll c and carotenoids. The width of the peaks around 
440 nm and around 675 nm varies from sample to sample, due to the change in 
accessory pigments present and the "package effect" [Morel and Bricaud 1981, Kirk 
1986]. The "package effect" and the accessory pigments/chlorophyll a ratio also cause 
the ratio ^(440) to chlorophyll concentration to be generally non-linear. 
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wavelength (nm) 

Figure 7d. Examples of Absorption Coefficients of Detritus. 

1.4. 1.4 Absorption Coefficient of Detritus . Figure 7d presents examples of the 
spectral absorption coefficient of detritus. Like gelbstoff, it is a strong absorber in the 
blue and exponentially decreases with wavelength, and can also be expressed as: 

a d = aj(440)e' S ^ X ' 4AO \ (10) 

The spectral slope S d has been reported in the range 0.009 to 0.011 nm' 1 [Roesler et 
al. 1989]. 

It has been found that the absorption coefficients for fulvic acid and 
detritus can be combined. The sum can also be expressed in a form similar to Eqs. 9 
and 10, but with an average slope at 440 nm of 0.011 nm' 1 [Carder et al. 1991]. 
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1.4.2 Volume Scattering Function and Scattering Coefficient 

The scattering property of a substance is determined by its volume scattering 
function, /3(a), which describes the probability of input photons being scattered to a 
specific angle, a. The scattering coefficient b is a measure of the sum of all the 
scattered photons (at all angles). 

1 .4.2.1 Volume Scattering Function. /3(a) . This property is defined as the scattered 
radiant intensity in the scattering angle a per unit scattering volume per incident 
irradiance. Figure 8a shows the VSF of sea water itself (/3 W ). As this is a molecular 
scattering mechanism, the scattering is symmetric between the forward and backward 
directions [Morel 1974]. 

Figure 8b shows the VSF of some natural waters [Petzold 1972]. Due to the 
presence of larger "soft," low-index particles, natural waters scatter photons strongly 
in the forward direction due to diffraction and refraction effects, and very weakly in 
the backward direction. 

1 .4.2.2 Scattering Coefficient, b . The integration of /3(a) over 4ir solid angle gives 
the total scattering coefficient. The scattering coefficient is divided into two parts: the 
forward scattering coefficient, b p and the back-scattering coefficient, b b . 
Mathematically they are 
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b f = 2itf o 2 p(a)sin(a)da, ^ 

b b = 2itJ*P(a)sin(<x)</a. 

1 

and b = b f + b h . b b is the major contributor to the upwelling radiance. 

The back-scattering coefficient of sea water can be divided into two 
components: back- scattering coefficient of pure sea water, b bw , and back-scattering 
coefficient of suspended particles, b bp . Everything except molecular scattering is 
included in particle scattering. 

Back-scattering coefficient of pure sea water, b bw . Many theoretical and 
experimental studies have been carried out [e.g. Morel 1974, Smith and Baker 1981] 
on the back-scattering of pure sea waters. Figure 9 shows the spectrum of the back- 
scattering coefficient of pure sea water [Smith and Baker 1981], which can be 
expressed as 

• (12) 

Due to the presence of ions and dissolved material, the scattering coefficient of 
pure sea water is -30% higher than that of pure fresh water [Morel 1974]. 

Back-scattering coefficient of particles, b kp . There are rare measurements of 
spectral b bp for ocean waters, although a few measurements at single wavelengths 
have been made (e.g. Petzold 1972), and Whitlock et al. [1981] reported b bp for river 
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wavelength (nm) 


Figure 9. Spectrum of the Back-scattering Coefficients of Pure Sea Water. 


samples. It is believed that b bp (k) can be expressed as [Morel and Prieur 1977, 
Bricaud et al. 1981, Smith and Baker 1981, Bricaud and Morel 1986, Gordon et al. 
1988, Sathyendranath et al. 1989a, Morel and Ahn 1990]: 

- vW 

with i) b in the range of 0 - 3 for different particles, while the data of Whitlock et al. 
[1981] imply an rj b of 1.7 for their river samples. 


400 y* 
* ; 


(13) 
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1.5 Relationships between AOPs and IOPs 


1.5.1 Sub-surface Irradiance Reflectance 


From Eq. 4, the sub-surface irradiance reflectance, /?(0"), is 


R( O') = 


EJL 0) 
Ejsn ' 


( 14 ) 


The major contribution to E d ( O') comes from the sun light and sky light. The 
contribution to £„(()'), however, is dominated by the following four components 
[Peacock 1992]: elastic scattering from molecules and particles in the water column 
(E u w (k)) and from bottom reflectance (E u b (X)), and inelastic scattering from gelbstoff 
fluorescence (E/(X)) and water Raman scattering (E u R (k)). Work by Gordon [1979] 
and Carder and Steward [1985] dealing with chlorophyll a fluorescence have been 
reported, but because the fluorescence efficiency of chlorophyll a varies by an order 
of magnitude and peak chlorophyll a fluorescence occurs in a narrow band centered 
around 685 nm, chlorophyll a fluorescence is not considered in this study. 

For optically deep, homogeneous water, the irradiance reflectance due to the 
elastic scattering of molecules and particles, R w (Qr), can be expressed as [Gordon 


et al., 1975]: 
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= £Lo G . 



(15) 


Recently, extensive Monte Carlo simulations have been performed for various sun 
angles and VSFs [Kirk 1984, 1991; Jerome et al. 1988, Gordon 1989b, Morel and 
Gentili 1991], and it is found that for bja < 0.25, generally, 


R w ( O') “ G(p/o)A, 

a 

where G(n d (0)), according to Kirk [1991], is 


(16) 


G(^0)) = Mj 





+M 2 , 


(17) 


M 1 and M 2 vary with VSF and range from 0.244 to 0.514 for M, and from 
0.303 to 0.326 for M 2 (for VSFs reported in Petzold [1972]). Actually, the above 
expression can be approximated as 


G(p^(0)) = A- , (18) 

P/0) 

for sun angles < 80° (0.68 < /T d (0) < 1.0) with G varying from 0.304 to 0.344, and 
averaged to 0.32 (±7%) for the above M, and M 2 values. Keeping G = 0.32 as a 
constant and putting the 7% error to the other terms, then Eq. 16 is 

/r(0‘) = 0.32D/0)— . (19) 

a 

This is a re-write of Eq. 16, but it more explicitly shows the variation of R w ( O') 


with solar zenith angle, as D d {Qi) = l//x d (0) ~ l/cos(/) (section 1.3.4). 
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For parts due to inelastic scattering and bottom reflectance, there are no simple 
expressions like Eq. 19. Their contributions to the remote-sensing reflectance will be 
discussed in Chapter 2. 


1.5.2 Remote-sensing Reflectance 


Based on the above assumption about the major components which make up the 
sub-surface upwelling light field, to first order (single scattering and quasi-single 
scattering [Gordon 1994]) the water-leaving radiance can be expressed as 


L = L : 


. b r f 1 R 
+ L„, + Li, + 


( 20 ) 


Breaking Eq. 6 into contributions from the various mechanisms listed in Eq. 20, 
we have 


Rrs + 


R f + 


Vs ' 


From Austin [1974], there is 


( 21 ) 


K • J ; i . (0 ' ) • <22) 

n, 

where t is the air-sea interface transmittance. Since E d (0 + ) = E d (0')/t, then from 
Eq. 6, 


R 


rs 


A L (O') A 

= i-RJ O') . 

nlEJLO) n 2 w 


(23) 
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Similar to Zaneveld [1982, 1994], through the radiative transfer equation, the 
exact solution for nadir remote-sensing reflectance from the elastic scattering part of 
the water column, R„ w ( O'), is 


where 


On 5 


c+k L -eb f 


( 24 ) 


f P( 

P* = -?■ . 05) 

J 2 L/B)cos(6)d( t } / 

is the light-averaged-backward-VSF (detailed in Appendix 1). 

The terms for bottom reflectance and inelastic scattering will be discussed in 
Chapter 2 and appendices. 


1.5.3 Diffuse-attenuation Coefficient of Downwelling Irradiance 


Through Monte Carlo simulation, Gordon [1989a] found that 

K/0) - IMDJMa+b,) , 
and 


K/av) « l.0SD/0)(fi+bJ , 


( 26 ) 


( 27 ) 


where /f^(0) and K d (av) are the diffuse attenuation coefficient for the downwelling 
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irradiance just beneath the surface, and an average of K d between the sub-surface E d 
and 10% of the sub-surface E d value, respectively. 
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CHAPTER 2 

MODEL OF THE REMOTE-SENSING REFLECTANCE 

2.1 Introduction 

One goal of ocean optics is to be able to interpret water color or remote-sensing 
reflectance in terms of the in-water constituents. From Eq. 21, we see that remote- 
sensing reflectance can be broken down into many terms. To interpret the measured 
remote-sensing reflectance on the left side of Eq. 21, each component on the right 
side of Eq. 21 needs to be expressed in terms of the optical properties of the water. 

2.2 Remote-sensing Reflectance of the Water Column 
(elastic scattering only), R rs w 

In Eq. 24, @ b = (3 bw + (3 bp , as /3 = 4- f3 p . (5 b is the light-field-weighted 

average of the backward VSF. Thus, (3 bw and (3 bp can not be simply expressed using a 
few parameters, since they depend on the light field. Extensive calculations for 
different a , VSFs, and light fields are necessary. 

Fory > 15° and a nadir viewing sensor, Gordon [1986] and Gordon et al. 

[1988] found through Monte Carlo simulation, that /?„ ,v (0') can be expressed as 
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C«n - EL 8 , 


h v 


a+b. i 

\ b J 


( 28 ) 


where g, = 0.0949, and g 2 = 0.0794. More recently, Morel and Gentili [1993] made 
Monte Carlo calculations for "case 1" waters ([chi a] < 3 mg/m 3 ) using one VSF 
shape for particles. For sun angles within 80° from zenith and for satellite viewing 
angles (0 a < 50°), they found 


On * g-> 

a 


( 29 ) 


where g averaged about 0.0936 (for 440 - 550 nm) and varied slightly with 
wavelength, sun angle and view angle. But, as pointed out by Mobley et al. [1993], g 
may vary with Monte Carlo computation models, and this variation can be 12% for 
high scattering waters and much larger in high absorbing waters [Mobley et al. 1993]. 
Also, it is not clear yet how g varies for high [chi a\ and "case 2" waters. 

In another way, with the definition Q = E u (0')/ L u (0) [Austin 1974] and Eq. 16, 
there is 


*;«o - 

Q a 


( 30 ) 


This expression introduces a new parameter, the Q factor. However, only a few 
measurements of Q exist. Q is extremely sensitive to sensor orientation, so its values 
have been reported from 3.2 to 12 [Gordon et al. 1980], For example, Austin [1979] 
takes Q to be about 4.7 and spectrally constant from 440 to 550 nm, while Kirk 
[1986] gives Q as -4.9, and Gordon et al. [1975, 1988] suggest a value of -3.4. In 
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some other studies, Q is somewhat arbitrarily chosen as a spectral constant [Carder 
and Steward 1985, Peacock et al. 1990]. However, recent measurements [C.O. Davis, 
unpublished] and Monte Carlo simulations for "case 1” waters [Morel and Gentili 
1993] show that Q increases with wavelength, which is an inverse trend compared to 
b bp [Carder et al. 1991]. 

Using the single and quasi-single scattering approximations of Gordon [1994], 
Lee et al. [1992, 1994a] found that for a wide variety of waters, Eq. 30 can be 
expanded as 



where Q m is the Q factor for molecular scattering. Q m can be estimated based on the 
molecular VSF shape and the illumination geometry [Lee et al. 1992, 1994a], 

Eq. 31 can be seen as a simplification of Eq. 24 for wide situations, with the 
first term in the bracket of the right side representing the backscattered photons due to 
molecular scattering, and the second term representing the backscattered photons due 
to particle scattering. 

In Eq. 31, both G and increase with sun angle [Gordon 1989b, Kirk 1991, 
Lee et al. 1992, 1994a], so we expect that /? rj (0) is less sensitive to the solar zenith 
angle than R(0 ). To simplify matters, we can take average values of G and Q m and 
allow any errors due to these simplifications to be embedded into the empirical terms 
X and Y. G is about 0.32 according to Gordon et al. [1975], Jerome et al. [1988], and 
Kirk [1991]. Q m is about 3.4 [Lee et al. 1992, 1994a]. Thus, our equation becomes 
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iC( O'A) 


0.32 

W*) 

a(X) 

3.4 



( 32 ) 


Variations due to VSF, solar zenith angle, and view angle are now all embedded in 
the two parameters X and Y. Since there is a different wavelength exponent between 
the scattering coefficients of molecules versus that of particles, errors generated by 
using average G and CL, values can not be completely compensated by the adjustment 
of the X and Y values. However, since b bw (\) is relatively small, the compensation can 
be close to 100% for turbid waters and larger than 90% for clear waters. 

Comparing Eqs. 29 and 32, it is found that if b bp of b b in Eq. 29 is expressed as 
X'(400/X)^, then X' * 3.42X + 0.0043/^(400). Since 0.0043^(400) is very small, 
Eq. 29 and Eq. 32 are mathematically almost the same. Thus, with two parameters to 
address b bp and/or Q, R rs w can be expressed by either Eq. 29 or Eq. 32. Using Eq. 

32, 3.42X can be used as an estimate of /> 6p (400). Also, since remote sensing 
measures part of the backscattered photons, Eq. 32 is more consistent with Eq. 24, 
the absolute solution. 


Actually, Eq. 29 could be derived from Eq. 24 if c + k L - e b f = 3 a, and = 
bj 3.55. However, both assumptions depend on the radiance distribution in the 
downwelling ( L d (Q,<p )) and upwelling (1.(0,?)) hemispheres. If the downwelling light 
is totally diffused, then J b = bjr. But this light field is not totally diffused in 
general, so (3 b ^ bjr. 

Eq. 32 pertains to optically deep water. In optically shallow water, scattering 


media and the backscattering signals are reduced due to the shorter water column. To 
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model these optically shallow waters, we consider that the subsurface E u w consists of 
two parts coming from two layers: one from the layer above the bottom, and one 
from the layer "below" the bottom. Then the subsurface E u w coming from the upper 
layer only can be obtained by reducing the optically deep expression by an amount 
equivalent to the contribution of the missing water column below H. Thus for shallow 
waters with depth H and a totally absorbing bottom, with t = 0.98, and n w = 1.34, 
R n w is approximated by 

RW _ O.nffyny J 400 > | y | 1 c 7.*<v r ->*| (33) 

" a [3.4 ( k ) } 

with z positive downward from the surface. 

If we define the quasi-diffuse attenuation coefficient as k = a+b h , then K d « 
D d K and K u * D u k [Gordon 1989a]. Since D u /D d « 2 [Gordon et al. 1975], we can 
rewrite Eq. 33 as 

R - . ml hi - 5u, n « 4 > 

where {D d } is the vertically averaged downwelling distribution function and {D d } is 
approximate by 1.080^0) [Gordon 1989a, Lee et al. 1994a]. 

2.3 Re mote- sen sing Reflectance of Bottom Reflectance, R rs b 


Assuming that the bottom is an extended Lambertian reflector with bottom 
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albedo p, then for a nadir-viewing sensor, R„ b can be approximated by [Lee et al. 

1992, 1994a] 

RZ (35) 

where k is the effective attenuation coefficient for the radiance from an extended 
Lambertian source. How k relates to the quasi-diffuse attenuation coefficient k is not 
well understood. Heuristically, it should be a value between c and k. It was found that 
k is from 1.36k to 1.62k for kH in the range of 0.5 - 2.0 based on the Monte Carlo 
simulations [Gordon 1989a] for an extended totally diffuse light source (detailed in 
Appendix 2). As an average in this work, k = 1.5k as is used by Marshall and Smith 
[1990]. Then Eq. 35 becomes [Lee et al. 1992, Lee et al. 1994a] 

- 0.17p«- W *' !w '. 0® 

This value depends not only on the optical properties of the water body, but 
also on water depth and the bottom albedo [Gordon and Brown 1974]. In the 
modeling work, the water depth was based on the Provisional Chart for the Gulf 
Coast (#1003), and the bottom albedo was based on measurements of bottom samples 
from the region with near-shore values of p = 0.1 to 0.2 (used for station TA01), and 
offshore values of p = 0.4 to 0.5 (used for stations TA02 and TA03). Figure 10 
shows examples of those albedo spectra. The quasi-diffuse attenuation coefficient k is 
assumed to be equal to total absorption a as b h < < a for most of the world ocean 
[Morel and Prieur 1977], 
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wavelength (nm) 

Figure 10. Examples of Bottom Albedo Spectra. 


For sensors not viewing the nadir, k in Eq. 35 needs to be adjusted to 
k! cos(9J. In our experiments, Q w is generally within 20°, i.e. 0.93 < cos(9J < 
1.0, so &/cos(0J * 1.5k could still be used. 


2.4 Remote-sensing Reflectance of Gelbstoff Fluorescence, R r f and 
Water Raman Scattering, Rj* 

In general, these terms are due to inelastic scattering. By considering (the 
volume scattering function for inelastic scattering) isotropic, remote-sensing 
reflectance due to gelbstoff fluorescence and water Raman scattering can be 
approximated as [Lee et al. 1992, 1994a] (detailed in Appendix 3) 
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, r K a(X x )E/0-,X x ) 

Ri(X) * 0.072 x](X x )- * * 


,+^r 


X [2fl(X)+a(X Jt )]£ 4i (0-,X) A 


-dX x , 


(37) 


and 


fl*(X) « 0.072 


b r (X x )EJ 0~,X x ) 
[2a(X)+a(X x )]EjL0-,X)‘ 


(38) 


So, combining Eqs. 34, 36, 37 and 38, after calculating R„*, Rj and R n b , only 
X and Y remain as unknowns. By matching the modeled R rs w and the residual of R n - 
Rj - Rj - R rs b , X and Y can be derived using the predictor-corrector approach to 
modeling as in Carder and Steward [1985]. 
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CHAPTER 3 
DATA AND METHODS 


3.1 Introduction 

From 1989 to 1993, measurements of optical properties for a variety of waters 
were taken, including sites in the North Atlantic, Monterey Bay, the West Florida 
Shelf, the Gulf of Mexico and the mouth of the Mississippi River. In these waters 
[chi <z] ranged from 0.07 - 50 mg/m 3 , and a g (440) ranged from 0.005 - 0.5 m"‘. 

Table 1 summarizes the cruise information, while Figure 2a shows the locations of 
data collection. For each station, hyperspectral R rs and particle and pigment absorption 
coefficients ( a p and aj of surface water samples were measured. For the 1993 Gulf of 
Mexico stations, a long-path (50 cm or 100 cm) [Peacock et al. 1994] 
spectrophotometer was used to measure a g . 

In-water optical properties ( K d and /?„(0» were calculated based on the profiles 
of downwelling irradiance and upwelling radiance, which were measured by Dr. J. L. 
Mueller (San Diego State University) or Dr. C. O. Davis (Jet Prop. Lab) using a 
Biospherical MER sensor (model 1048 A). 
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3.2 Remote-sensing Reflectance, R rs 


Hyperspectral R n was measured by the method developed by Carder and 
Steward [1985], using a Spectron Engineering spectroradiometer (Spectron Model 
SE-590) with 253 spectral channels covering the wavelength range from 370 - 1100 
nm. With this instrument, the upwelling radiance above the surface (L U (Q + ,Q a ,tp)) was 
directly measured, with 0 a < 30° and about 90° from the solar plane (see Figure 5 
for the geometry). The downwelling sky radiance (L sky ) was also directly measured in 
the same plane as L u ( 0 + ) but from a direction reciprocal to L„(0 + ). Downwelling 
irradiance was derived by measuring the radiance (L G ) reflected from a standard 
diffuse reflector (Spectralon). With these measurements, R rs was derived through 


L(0*)-rL. 

R = — !!- — ^R g - A . 


7 lL r 


(39) 


The upwelling radiance is 

i.(0-) - (40) 

where r is the Fresnel reflectance of the water surface, and R a is the reflectance of 
the diffuse reflector. Figure 1 1 shows the values of sea water r for the horizontal and 
vertical electric components as well as for the unpolarized light, r ~ 0.018 (a 20° 
view angle) is usually used when there is a vertical polarizer in front of the sensor, 
and A is the offset for sun glint due to waves and surface foam, which is derived by 
forcing /?„(750) = 0 (usually for offshore waters). 
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Figure 11. Fresnel Reflectance of Sea Water ( n w = 1.34). 


3.3 Absorption Coefficient of Particles and Pigments, a p and a ^ 


The method described by Mitchell and Kiefer [1988] was used to measure a p , 
and the method developed by Kishino et al. [1985] and modified by Roesler et al. 
[1989] was used to measure a d and to derive a 4> . Briefly, water samples were filtered 
onto Whatman GF/F filter pads immediately after collection. Another pad wetted with 
filtered sea water served as a blank. The transmission spectra from 380 nm to 800 nm 
of these pads (T p and T b ) were measured by the Spectron. The optical geometry was 
designed to illuminate the pad with diffuse light from a Lambertian diffuser added 
between the light source and the filter pad. This geometry is very similar to that of 
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Table 2. Recent Results for "Beta Factor" of GF/F Filter Pad. 



Bricaud and Stramski [1990], who illuminated a pad in front of a diffusing window 
which was adjacent to an end-on photomultiplier tube. The volume filtered was about 
500 - 1500 ml (depending on the water clarity). 

Using the measured transmission spectra, the optical density of the sample is 

OD p = \o%— , < 41 > 

p T 

p 

and the absorption coefficient of the sample is 


OD b 

a = 2.3 p - 

P 


where B is the ratio of the volume filtered to the effective surface area of sample on 


the pad. is the so called "beta factor" [Mitchell and Kiefer 1988], introduced to 
describe the optical-path-length elongation due to the filter pad. Table 2 and Figure 12 
provide recent results of ^ for GF/F filter pad. Notice the significant differences of 
when OD p is around or less than 0.1. In our calculation of a p , (3^ from Bricaud 
and Stramski [1990] (their Eq. 2) was used. In the calculation of the absorption 
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Figure 12. Recent Results for "Beta Factor" of GF/F Filter Pad. 

coefficient, large-particle scattering was corrected by assuming a p (780) = 0 from 
the a p curve. 

After this measurement, the sample pad was soaked in hot methanol [Kishino et 
al. 1985, Roesler et al. 1989] for about 15 minutes to remove pigments, and its 
optical density was again determined via Eq. 41, and the detrital absorption 
coefficient a d was obtained by Eq. 42. 

The difference between the particle and detrital absorption coefficients provided 
the absorption coefficient of phytoplankton pigments, a 



( 43 ) 
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3.4 Absorption Coefficient of Gelbstoff, a g 

For stations before 1993, the absorption coefficient of gelbstoff (a g ) was not 
explicitly measured for visible wavelengths. However, from Eq. 27, the total 
absorption coefficient could be derived from K/av), and a w and a p were available, so 
a g could be estimated from K d {av) values, using the expression a g * 

AT rf (av)/(l .080^(0)) - a„ - a p , as b b < < a. At the 1993 GOMEX and COLOR stations, 
a g of water samples were measured using 50-cm or 100-cm path-length instruments, 
respectively, after filtering the sample through 0.2 /im pore-diameter Gelman Supor- 
200 filters [Peacock et al. 1994], 

3.5 Measurements of In-water Optical Properties, K d and R rs (0) 

A Biospherical Instrument MER (model 1048A) was used to determine the 
vertical structure of the water column. The depth profiles of down welling irradiance 
( E/z )) and upwelling radiance (L/z)) were determined for each cast. 

With these measurements, K d {av) and sub-surface downwelling irradiance 
(E/O')) and upwelling radiance (L u ((f)) were derived. The process is similar to that of 
Smith and Baker [1981]. Briefly, 

-K/av)z 


E/z) = E/0')e 


(44) 
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then, 


In [Ejtz)] = ME/0-)]-Kjtav)z . (45) 

Due to surface-wave-focusing and ship shadow effects [Gordon 1985], the 
measured E d (z) and z include errors associated in the field, and we cannot simply use 
values from two depths to derive K/av) and E d { O'). In order to correct these errors, 
linear regression between ln[£^(z)] and z were performed for the surface layer. The 
regression results gave K d (av) and ln[£ d (0)]. The same process was applied to L u (z) to 
get L u { 0 ). 
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CHAPTER 4 

RESULTS OF R„ MODELING 
4. 1 Validation of R n Model 

In previous studies [Carder and Steward 1985, Peacock et al. 1990, Lee et al. 
1992, Lee et al. 1994a], qualifying words such as "excellent," "very good” or "good" 
were used to describe how well the modeled curve fit the measured curve, and they 
were usually justified visually. There was no quantitative indication regarding the 
difference between the measured and modeled curves. Here, a modified average 
percentage difference ( a.p.d .) is used to quantify and validate the model results of R rs . 
It is defined as 

a ' P ' ' AVG™(R™)+AVG™aiZ“) 

where AVG^ means the average value in the wavelength range from to X 2 . 

The cut-off between 660 and 750 nm is because there is no term included to 
model the chlorophyll a fluorescence in the measured signal. The 750 - 830 nm band 
is important for turbid waters. 

In the forward process, a g was measured or estimated, a p was measured, H for 
shallow waters came from published charts, and the predictor-corrector method 
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Table 3. Station Chosen for the R rs Model Demonstration. 


Station 

Latitude 

Longitude 

Time/Date 

Bottom 

depth 

Mod. 

depth 

TA01 

27°27’ N 

82°55’ W 

+ 10.5/3-4-90 

14m 

13.7m 

TA02 

27°20’ N 

83°03’ W 

+ 13. 0/3-4-90 

25m 

25m 

TA03 

27°12’ N 

83°ir W 

+ 14.9/3-4-90 

35 m 

36m 

GO08 

28°48’ N 

91°30’ W 

*08.5/4-12-93 

- 

- 

GO10 

28°15’ N 

91°30’ W 

*14.0/4-12-93 

- 

- 

G027 

29°32’ N 

85°47’ W 

‘09.2/4-19-93 

- 

- 

C012 

28°52’ N 

89°33’ W 

+ + 10. 8/6-5-93 

- 

- 

COM 

28°48’ N 

90 < X)2’ W 

+ + 16. 1/6-5-93 

- 

- 

C019 

27°34’ N 

83°20’ W 

+ + 09. 5/6-8-93 

33m 

35m 


Note: " - indicates that the water is optically deep. + : East standard time, 
*: Central daylight time, + + : East daylight time. 


[Carder and Steward 1985] was used to obtain the values of X and Y. In this process, 
it was sometimes necessary to apply a factor/ to the measured a p , such that the 
absorption coefficient of particles used in the model is /* a p . The purpose of this /is 
to correct for possible errors due to patchiness, and/or the "beta factor" (see Figure 
12). It is necessary to keep in mind that the process here was focused on testing the 
expressions developed in Chapter 2 using the measured components, so only minor 
adjustments were applied to the measured components before a small (< 0.05) a.p.d. 
was obtained. Most of the parameter derivation was concentrated on X and Y (and p 
for shallow waters), and the parameter / was not varied much (usually 0.9 to 1.1) 
from 1.0 in reducing the a.p.d.. 
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For the calculation of Rj, n, s, )y and a came from the measurements of Hawes 
et al. [1992] and Hawes [1992], while for R„* came from Collins et al. [1984]. 

4.2 Model Results and Discussion 

Measured R n curves for 9 stations were selected to demonstrate the model 
developed in Chapter 2. These 9 stations covered a wide range of water types: 1) the 
West Florida Shelf with shallow, gelbstoff-rich coastal waters; and 2) Gulf of Mexico 
waters with phytoplankton blooms in the Mississippi River plume (S > 17%o). 

Table 3 provides the station locations as well as the measured and modeled water 
depths for the shallow stations. Figure 2b shows the locations of these stations in the 
Gulf of Mexico. As examples, Figures 13a - 13d show the detailed model components 
for R rs , and Figure 14 shows the results for the chosen stations. Table 4 lists the 
model parameters X, Y, a g ( 440), and 0,(440) with a.p.d. and the measured values of 
0,(440) and [chi 0 ] for each station. Table 5 details the fractional contributions that 
R n w , Rj, Rj and R n b make to the measured R rs at 440 nm and 550 nm. 

It can be seen from Figures 13a - 13d that excellent fits were achieved between 
the measured and modeled R n curves for all of the selected stations except for the 
spectral region around 685 nm where chlorophyll 0 fluorescence is present in the field 
data. This excellence can also be seen from the a.p.d. values (Table 4) for each 
station with the average a.p.d. being 3.1%, which is within the measurement 
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Figure 13a. Detailed R n Model for Station TA01. 





wavelength (nm) 


Figure 13b. Detailed Model for Station G027. 
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accuracy. For such a small percentage difference, we can reasonably say the model 
developed here works very well for these waters. 

For the stations chosen, the a p that was required by the model was in general 
within 10% of the measured a p except near the Mississippi River at station C012 
(20%), with very high [chi a] (38.6 mg/m3). The average difference between the 
measured and required a p is 8.9% (9 stations) (7.5% when station C012 is excluded). 
The maximum 15% or 20% difference can perhaps be explained by the water 
patchiness and/or the accuracy involved in the method of a p measurement due to the 
"beta factor" which varies significantly between species and researchers [Bricaud and 
Stramski 1990, Mitchell 1990, Yentsch and Phinny 1992, Cleveland and Weidemann 
1993, Nelson and Robertson 1993] (see Table 2 and Figure 12). This effect may be 
especially important for station CO 12 which was near the Mississippi River mouth 
where the heavy load of sediments and minerals might cause additional uncertainty in 
the optical path length elongation. Also, the influence of horizontal and vertical 
structure of the waters increases for mesotrophic - eutrophic waters, so patchiness can 
affect accuracies in the more hypertrophic waters. Finally, the low signal obtained for 
the upwelling radiance measurements at station CO 12 made the R rs calculation 
sensitive to corrections for reflected skylight. 

It can be seen from Table 4 that the ratio of 0,(440) to o p (440) among these 
waters was highly variable, with a range from 0.3 to 3.0, and the X value does not 
co-vary with the pigment concentration [chi a] for the waters studied. This illustrates 
that the model works well over a wide range of conditions, and also suggests why the 
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Figure 14. Measured vs. Modeled R n for the Selected Stations. 
Solid: measured, dashed: modeled. 


power-law pigment algorithm does not work well for coastal waters because of the 
lack of co-variation of all optical components with [chi a]. The highest X value, 
0.0087 m 'sr 1 , was at the shallow, mesotrophic station TA01, suggesting a large 
scattering effect due to detritus and suspended sediments. Brisk northwesterly winds 
suspended sediments in the shoal regions to the east and north of the station, and they 
likely were transported by the ebb tidal currents from Tampa Bay to the study site 

[Carder et al. 1993]. 

The Y values were within the range from 0 to 2.4 for the waters reported here, 
which might be interpreted as being partially due to b bp . For b bp , as mentioned m 
section 1. 4.2.2, the wavelength exponent (tj„) is in the expected range 0 - 3.0 for a 
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Table 4. Model Parameters for Stations in Table 3. 


Station 

X 

Y 


a P> 

a.p.d. 

men 

°pl 

[chi a ] 

TA01 

.0087 

1.2 

.082 

.041 

.037 

.045 

1.05 

TA02 

.0022 

2.4 

.042 

.033 

.025 

.035 

.61 

TA03 

.0012 

2.3 

.034 

.027 

.029 

.026 

.70 

GO08 

.0065 

.24 

.31 

.28 

.033 

.295 

5.27 

GO10 

.0010 

1.8 

.059 

.023 

.023 

.021 

.12 

G027 

.0014 

1.9 

.078 

.029 

.021 

.034 

.20 

C012 

.0030 

0 

.42 

1.34 

.044 

1.12 

38.58 

C014 

.0068 

0 

.38 

1.21 

.031 

1.16 

20.26 

C019 

.0007 

2.0 

.023 

.021 

.037 

.023 

.22 


Note: a gl = a g ( 440), a pI = g/440). 


range of particles (e.g. river sample [Whitlock et al. 1981], bacteria [Morel and Ahn 
1990], phytoplankton cells [Bricaud and Morel 1986], and coccoliths [Gordon et al. 
1988]). For Y = 2.4, if the exponent for Q is 1.0, then r\ b = 1.4, which is within the 
0 - 3.0 range reported elsewhere, also is within the range of 0 - 2 as used by 
Sathyendranath et al. [1989a]. Generally Y values were found to be low for turbid 
water, and high for clear water. The value 2.4 for stations TA02 and TA03 of the 
West Florida Shelf seems a little bit high for those waters. This might be due to 
errors in a g , as a g was derived from K d for those stations. For the world ocean, the 
range of Y needs further study. 

At the optically shallow stations (TA01, TA02, TA03, and C019), the model- 
required depths were within about 10% of the chart depths without consideration of 
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any tidal influence (typically <0.5 m). This demonstrates that Eq. 36 works well for 
shallow waters, and it indicates a potential to use this model to survey (e.g. by 
aircraft overflights) dramatic changes in shelf bathymetry that can occur as a result of 
major storms. For station C019, an albedo value of about 0.1 was used in the model, 
suggesting the bottom might contain more heavy minerals or grass at that site. Direct 
bottom albedo measurements are lacking at individual stations and are needed for a 
wide variety of bottom types. 

For stations TA01, TA02, TA03, GO08, C012 and C014, the general 
agreement between the modeled and measured R n values is very good ( — 3.5% 
a.p.d.), with small differences around 580 nm, where the measured R n > modeled 
R rs . Other than modeling error, there are at least three possible reasons for this: a) 
bottom albedo uncertainty, b) phycoerythrin fluorescence [Yentsch and Yentsch 1979, 
Yentsch and Phinney 1985], and c) water absorption coefficient uncertainty [Smith 
and Baker 1981, Tam and Patel 1979]. 

A spectrally constant bottom albedo was used for the shallow stations. Soil 
reflectance [Tucker and Miller 1977] and earlier measurements of bottom albedo 
(Figure 10) did display some spectral dependence. The amount of change spectral 
albedo could induce would not provide the sharp spectral increase and then decrease 
with wavelength in R n required for the measured and modeled R rs curves to converge, 
however. There also was no bottom contribution to R rs at GO08, CO 12 and COM. 
More realistic explanations include the lack of a term in the model for phycoerythrin 
fluorescence, the differences between the water absorption coefficients in this 
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Table 5. Optical Component Contributions to R„ of Stations in Table 3. 


Station 



Rn/Rrs 

RJ440)/RJ550) 

440 

550 

440 

550 

440 

550 

mea. R n 

corr. R n 

TA 01 

.95 

.92 

.02 

.01 

.02 

.10 

.87 

.91 

TA02 

.91 

.82 

.05 

.06 

.05 

.11 

1.48 

1.77 

TA03 

.87 

.83 

.10 

.12 

.02 

.04 

1.87 

1.96 

GO08 

.97 

1.00 

.03 

.02 

- 

- 

.49 

.48 

GO10 

.92 

.94 

.08 

.07 

- 

- 

1.64 

1.61 

G027 

.91 

.96 

b 

OO 

.07 

- 

- 

1.34 

1.26 

C012 

.92 

1.03 

.09 

.03 

- 

- 

.31 

.28 

C014 

.99 

1.02 

.03 

.01 

- 

- 

.32 

.31 

C019 

.92 

.96 

.06 

.07 

.01 

.01 

2.46 

2.37 


Note: R* = Rj + R„ r ; corr. R n = mea. R rs - R* - R n b . 


spectral region as reported by Smith and Baker [1981] and by Tam and Patel [1979], 
and the accuracy of the "beta factor" for low absorbing portions of c^(X) for detritus 
-rich stations. Further study is required in order to resolve this issue. The differences 
between the measured and modeled R rs curves around 685 nm, on the other hand, are 
expected due to the fact that no term is included in the model to describe the 
chlorophyll a fluorescence [Carder and Steward 1985], 


4.3 Contributions by Rj and R„ r 


Model results at station TA03 suggest relatively higher gelbstoff fluorescence 
and water Raman scattering influences, since a higher 77 (1.5%) [Hawes et al. 1992] 
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was used. This value is — 3 times greater than the value suggested by Spitzer and 
Dirks [1985] for terrigenous gelbstoff. If we exclude this station, 90% or more of the 
water-leaving radiance is accounted for by the sum of the elastic scattering from 
molecules, particles and the bottom, which leaves about 10% or less of the measured 
/?„ for gelbstoff fluorescence and water Raman scattering. This is consistent with the 
reports of Marshall and Smith [1990] and Stavn [1990], as water Raman scattering 
makes more of a contribution when the water is clear. 

It is interesting that the ratio R rj (440)//? rj (550) did not vary widely due to 
inelastic scattering (see Table 5). Among stations without bottom influence, 
differences in the ratio were within - 10%, which suggests the spectral radiance ratio 
algorithm is effective for most deep waters without consideration of gelbstoff 
fluorescence and water Raman scattering. However, it is obvious that as the bottom 
influence increases, the usefulness of the power-law algorithm decreases. Also, the 
power-law algorithm can not distinguish between the absorption of gelbstoff and that 
of pigments. Note that J?„ w (490)/J? ff (490) values as low as 0.77 were determined (not 
explicitly shown), suggesting that great care must be taken when interpreting remote- 
sensing curves for the intermediate wavelengths at shallow coastal stations. 
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CHAPTER 5 

CONCLUSIONS OF PART A 


5.1 Contributions to the water-leaving radiance spectra for a variety of waters 
can be attributed to elastic scattering by water molecules, suspended particles, and 
bottom reflectance, and to inelastic scattering by water Raman scattering and gelbstoff 
fluorescence. Inelastic scattering by pigments was not considered. For optically deep 
water, remote-sensing reflectance of the water column part (elastic scattering only), 
R rs w (X), can be mathematically simulated as follows 




0.17 


<* w (\)+a g (\)+a p (\) 

3.4 



(47) 


where b bw (k) is known, X and Y are spectral constants, and Y was less than 2.4 for the 
waters considered. For optically shallow waters, the expression for bottom-reflectance 
contribution, 


R* « 0.l7pe' (lDj * l5)aH , (48) 

works well for the shallow waters considered. Together, the water-column term and 
the bottom-reflectance term accounted for about 90% or more of the total remote- 
sensing reflectance. 
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5.2 Close agreement between modeled and measured R n was achieved for all 
selected stations when all of the scattering mechanisms mentioned (both elastic and 
inelastic) were included. The average a.p.d. is 3.1%. The ratio a i (440)/a /> (440) 
ranged from ~0.3 to -3.0, indicating the broad usefulness of the model. For 
contributions other than from the water column, as much as 23% of f?„(490) is 
attributable to water Raman scattering, gelbstoff fluorescence, and bottom reflectance 
for an optically shallow (25 m) station. For pigment algorithms based on the power 
law of spectral radiance ratio, most "error" comes from reflected bottom radiance for 
optically shallow coastal waters. 

The a p required by the model was generally within 10% of the measured a p with 
an average difference of 8.9% (9 stations) (7.5% when station C012 is excluded). 

This suggests a potential to remotely measure the pigment and gelbstoff absorption 
coefficients, although derivation of [chi a] will depend upon knowledge of the 
chlorophyll-specific absorption coefficient for a region. 

5.3 The model-required bottom depths for the optically shallow waters are 
within 10% of the chart depths, suggesting its possible use to remotely measure 
bottom depth for the shelf waters. 

5.4 The contribution of Rj and R rs R were generally within 10% or less of the 
total R„, covering the whole range from 400 nm to 600 nm. If these two terms were 
omitted, this 10% or less difference can be roughly compensated by a small increase 
in X value in the modeling. Also, Rj and R a * do not significantly affect the 
/?„(440)//?„(550) ratio; thus, the power-law pigment algorithm can be used without the 
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correction of gelbstoff fluorescence and water Raman scattering with little error if the 
absorption and scattering properties co-vary with pigment concentration as happened 
for "case 1" waters. These imply that Rj and R r * could be omitted in the R n 
modeling to simplify the process. 
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PART B 


THE INVERSE PROBLEM AND APPLICATIONS 
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CHAPTER 6 

THE INVERSE PROBLEM 


6.1 Introduction 

As discussed previously, the ultimate goal of remote sensing is to derive the in- 
water constituents through remotely measured signals. Empirical and semi-analytical 
approaches have been discussed for derivation of pigment concentration [Clark 1981, 
Gordon and Morel 1983, Carder et al. 1991, Carder et al. 1994], and diffuse 
attenuation coefficient [Austin and Petzold 1981, Gordon and Morel 1983]. For the 
analytical approach, an inverse of the process described in Chapter 2, not much 
improvement has been achieved in the past studies. One of the difficulties is that a ’ 
or aj change drastically from region to region, and it is not easy to accurately 
express these changes by a few parameters such as for a g or a d . 

For R„ of deep waters at N wavelengths (A„ \ 2 , .... \ N ), ignoring Rj and R a * 
as discussed in section 4.3, there are N equations 
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in which a dg can be expressed as [Roesler et al. 1989, Carder et al. 1991] 


(49) 


-*») - V 440) e - 5 *"-» . (50) 

There are at least N+4 unknowns for the N equations (N for a<,(X), 2 for a dg (k ) 
(a dg ( 440 ) and 5^) and 2 for particle scattering (X and F)). If only R rs is available, 
there will be no certain solution for the above equations unless we dramatically reduce 
the unknowns regarding <2*(X). 

Bidigare et al. [1990] pointed out that a* can be re-constructed by knowing 
concentrations and the specific absorption coefficients for each pigment, but these 
cannot be known from remotely sensed data. Hoepffner and Sathyendranath [1991] 
suggested that can be modeled by the sum of 1 1 Gaussian bands. For these 1 1 
Gaussian bands, their center wavelengths and half band-widths would vary from 
phytoplankton species to species. Even if the center wavelengths and half band-widths 
can be fixed, we still need 1 1 parameters to simulate a*. Methods are also suggested 
to use average specific absorption coefficient [Morel 1980, Sathyendranath and Platt 
1988, Carder et al. 1991] or average absorption curves [Roesler and Perry 1994]. 
With these approaches, only one unknown (the pigment concentration or a scale 
factor) is needed to model a*. Thus, if N is equal to or greater than 5, theoretically 
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the series of N equations could be solved and the unknowns related to the absorption 
and scattering could be derived. But, due to the "package effect" and changing 
environments, it is well known that the curves vary widely from sample to sample 
[Morel 1980, Bricaud and Stramski 1990, Bidigare et al. 1990, Hoepffner and 
Sathyendranath 1991, Carder et al. 1994], No single shape or value for aj can be 
used globally. So, for the inverse problem in remote sensing, simpler expressions 
with adequate accuracy for a * would be very useful. The following section will 
discuss the possible simple methods to simulate a^X) with the consideration of the 
change of a ^ shape with water sample. 

6.2 Simulation of a^X) 

Generally, there are two ways to simulate a*(\): one is by the combination of 
mathematical functions [Hoepffner and Sathyendranath 1991, Lee et al. 1994b], which 
use a few parameters to simulate the whole spectrum; another is to empirically relate 
a* of each wavelength to a specific value such as total pigment concentration [Morel 
1980], or [chi a] [Carder et al. 1991] or at one wavelength [Carder et al. 1994, 
Roesler and Perry 1994], The first method is simple and has more power to adjust the 
a* curve. This method, however, creates a smooth a * curve, and sacrifices the finesse 
of the a* curve containing the pigment composition information when the number of 
simulation parameters is limited. The second method shows the averaged finesse of 
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the a ^ curve, and has less potential to adjust the curve as generally it is controlled 
by one parameter. 


6.2.1 Mathematical functions 

By analyzing surface d^(X) data collected from the Gulf of Mexico in April, 
1993, which covered a [chi a] range from 0.07 to 40 mg/m 3 , an expression for a/X) 
was suggested by Lee et al. [1994b], which is a combination of 3 simple functions 
involving 6 parameters. Among the 6 parameters, 2 parameters vary only slightly for 
different waters, and only 2 parameters have strong effects on the whole a * curve. 
There were a few a * curves from deep water samples which were greatly different 
from those of the surface samples, and could not be well simulated by the suggested 
simple expressions. But this is not important for remote sensing as 90% of the 
observed photons derive from the top attenuation depth ( 1 / K d ) [Gordon and McCluney 
1975]. For species recognition, the methods of Bidigare [1990] and Hoepffner and 
Sathyendranath [1991] might be better. 

For the wavelength range of 400 nm < X < 700 nm, the simple mathematical 


simulation for < 2 *(X) is: 
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400 s X s 570, ^(A.) = a^e 


( 51 ) 


a.(656) - a. (570) .... 

570 < X < 656, aJX) = a A (570) + -* * (X - 570), < 52 > 

♦ ♦ 656 - 570 


and 


(x - xj? 

656 5 1 s 700, a/X) = a^e 2o * 


(53) 


The wavelength range for Eq. 51 was 400 - 590 nm in Lee et al. [1994b], but it 
is better for R rs modeling if the range is adjusted to 400 - 570 nm. 

With Eqs. 51 - 53, a 9 curves can be simulated. Figures 15a - 15d show 
examples of the simulated versus measured a For the a * samples, the normalized 
root-mean-square (rms) error is in the range of 5 - 20% with an average of 11%. 

Most of the variation occurred around 570 nm for clear water stations when the 
measured a * values were low. 

In Eqs. 51 - 53, there are 6 parameters, a <t>1 , F, a <j)2 , X, and o 2 . Parameter F 

describes the width of the a ^ curve from 400 nm to ~560 nm, 100 -I- X ; is the 
wavelength of the blue peak, X 2 is the wavelength of the red peak and 2.355a 2 
determines the half band-width around the red peak. For the samples studied [Lee et 
al. 1994b], F varies from 1.6 to 4.2, X, varies from 338 - 342 nm with 80% at 340 
nm, X 2 varies from 672 to 675 nm with most at 674 nm, and 2.355a 2 ranged from 21 
to 34 nm. a <t>1 varies from 0.01 to 0.83 m' 1 while a^/a^, varies from 0.21 to 0.85. So, 
for the 6 parameters, and X 2 are almost fixed, and a^ 2 and <x 2 only affect a small 







Figure 15c. Measured vs. Simulated a 6 for Station G015. 



wavelength (nm) 

Figure 15d. Measured vs. Simulated a* for Station CQ15. 
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Figure 16a. Parameter Fof Simulation vs. <2^(440). 
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Figure 16b. Parameter a 2 of Simulation vs. <2^(440). 
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Figure 16c. Ratio a*(675)/a*(440) vs. o*(440). 

Solid curve: Eq. 54; filled circle: Eq. 55. 


region of the a * curves. Thus, only 2 parameters, a H and F, are important to the a* 
curve. 

F, <j 2 and are indicators of the "package effect". The greater the 
package effect, the smaller the parameter F, the bigger the band-width 2. 355a, and 
larger the ratio aja^. That means, for in vitro phytoplankton pigment absorption 
coefficients (i.e. no "package effect"), from the above results, the "fatness" factor F 
will be close to 4.2, around 0.2, and the half band-width (2.355a,) around the 
red peak will be close to 21 nm (a similar value as reported by Hoepffner and 
Sathyendranath [1991]). 

Unlike Lee et al. [1994b] who related the parameters to measured R n , the 
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parameters can be related to a*(440) after nonlinear regression analysis. Figures 16a - 
16c show how those parameters relate to the measured a*(440). It is found that: 
a d a *i = 0.86 + 0.161n(^), normalized rms error: 17.2%, 

F = 2. 89exp[-0.505tanh[0.561n(a^/0. 043)]], normalized rms error: 12.4%, 

0 2 = 14.17 + 0.91n(a^), normalized rms error: 5.6%. 

In this way, the a * curve can still be estimated if the measured R n is contaminated by 
bottom reflectance. 


6.2.2 Empirical Relationship 

For the measured surface a ^ of cruises GOMEX and COLOR, Carder et al. 
[1994] found that for the SeaWiFS channels, a^(X)/^(675) can be expressed by a 
hyperbolic tangent function, i.e. 

a^X) = a 0 (X)e a ‘ (i) ^ (XMa * (675)/ai(i))1 ^(675) , < 54 ) 

where the parameters a 0 , a,, a 2 , and a 3 are empirically determined for each SeaWiFS 
wavelength [Carder et al. 1994]. By the above expression, the value of fl 4 ,(X)/a <> (675) 
approaches an asymptote at very high or very low values of <3^(675) (see Figure 16c). 

For 0^(440) in the range of 0.01 to 1.0 m 1 ([chi a\ equivalent is in the range of 
0.07 to 50 mg/m 3 ), a simplified expression is 

V*) 


V440) 


= a 0 (X) +a 1 (X)ln[a <J( (440)] , 


(55) 
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with parameters a 0 and a t empirically derived for each wavelength and presented in 
Table 6. 


Table 6. Parameters for the Empirical a*(X) Simulation. 


X 

a 0 


1 x 

a 0 

0 / 

r 2 

390 

.5813 

.0235 

.057 

560 

.3433 

.0659 

.784 

400 

.6843 

.0205 

.081 

570 

.2950 

.0600 

.806 

410 

.7782 

.0129 

.058 

580 

.2784 

.0581 

.834 

420 

.8637 

.0064 

.032 

590 

.2595 

.0540 

.848 

430 

.9603 

.0017 

.005 

600 

.2389 

.0495 

.845 

440 

1.0 

0 

— 

610 

.2745 

.0578 

.875 

450 

.9634 

.0060 

.113 

620 

.3197 

.0674 

.892 

460 

.9311 

.0109 

.220 

630 

.3421 

.0718 

oo 

\o 

OO 

470 

.8697 

.0157 

.235 

640 

.3331 

.0685 

.893 

480 

.7890 

.0152 

.179 

650 

.3502 

.0713 

.891 

490 

.7558 

.0256 

.356 

660 

.5610 

.1128 

.884 

500 

.7333 

.0559 

.710 

670 

.8435 ; 

.1595 

.893 

510 

.6911 

.0865 

.815 

680 

.7485 

.1388 

.886 

520 

.6327 

.0981 

.836 

690 

.3890 

.0812 

OO 

O 

530 

.5681 

.0969 

.823 

700 

.1360 

.0317 

.751 

540 

.5046 

.0900 

.805 

710 

.0545 

.0128 

.645 

550 

.4262 

.0781 

.779 

720 

.0250 

.0054 

.531 


Using these expressions (Eqs. 54 or 55), the number of unknowns related to 
a^(\) is reduced to 1 (<z 0 (675) or a 0 (44O)). However, with 4 or 2 parameters for each 
wavelength, the shape will no longer be the same for different waters, and the 
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changes of a*(X) due to "package effect" or pigment composition are considered, at 
least to the first order. 

Due to the limited data sets, it is improper to claim that the values in Table 6 
are final and universally usable. However, these values work very well for the waters 
in this study. Improved empirical results are anticipated with increased data sets and 
better measurements. 

In the following sections, simulations of a^X) by mathematical functions (Eqs. 
51 - 53) and by empirical relationships (Eq. 55) are tested in the R rs inversion for 
waters from the Monterey Bay and Gulf of Mexico, where both measured R rs (X) and 
K d (X) are available. In the application chapter, efforts are concentrated on using the 
mathematical simulation, however. 

6.3 Derivation of the Absorption Coefficient from R n 

From the above a ^ simulations, it is found that 440) is the most important 
parameter for the a^X) curve. The other parameters can either be fixed or be 
estimated from <^(440). For instance, for the mathematical simulation method, X, and 
X, can be fixed at 340 nm and 674 nm, respectively. If the relationships of a^/a^, a 2 
and F versus a ^ are used in the R rs inversion, the number of total unknowns is 
reduced to 5 for the N equations: a#, a dg { 440), S dg , X and Y. By minimizing the 
a.p.d. defined in section 4.1, it is possible to derive the 5 unknowns if N > 5. 
Before we do this, however, ranges for the unknowns have to be set as there exist 
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realistic limits for them. It is possible that values outside of these limits may provide 
smaller a.p.d. values. 

For the exponent Y, there are no measurements available. Part of Y is rj b , and as 
discussed in section 1.4.2, rj b changes with particle size. Generally, it is assumed that 
Vb ~ 10 for open ocean and that n b ~ 0 for coastal waters [Gordon and Morel 1983, 
Morel 1988], although rj b can be as high as 1.7 for river samples [Whitlock et al. 
1981] and 3.0 for coccolithophorid blooms [Gordon et al. 1988], Due to the similar 
curvature of the b b and a dg spectra, however, the range for Y cannot be simply set as 
0 — Y ^ 3, because when the absorption is dominated by a dg (very common for 
coastal waters), the compensation between the a dg and b b parameters becomes strong. 
Therefore a narrow range for Y for each station must be specified. Previous model 
results yield a rough relationship similar to that of Carder et al. [1994]: Y * 0.86 + 
1.21n(x) with x = R ra (440)//? ra (490) . Thus, the range for Y is set as: 

0.9 * (0.86+ 1.21n( x )) < Y< 1.1 * (0.86+ 1.21n( x )), 
i.e., within 10% of the regression value and keeping T > 0. 

S dg , which depends on the relative abundance between detritus and gelbstoff, 
vanes from sample to sample [Roesler et al. 1989], By considering a detritus-to- 
gelbstoff ratio less than 1.0, the range for S dg is set as 
0.012 < S dg < 0.016, 

as S g varies from 0.011 nm' 1 to 0.019 nm' 1 for different materials [Carder et al. 1989, 
Hawes et al. 1992], and averaged about 0.014 nm 1 for ocean waters [Bricaud et al. 
1981]. 
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The ranges for a^, a dg and X are much easier to establish: 

> 0, a di (440) > 0, and X > 0. 

By minimizing the a.p.d. defined in Chapter 4, the 5 unknowns are derived for each 
measured R n curve. Since there were R rs values at —180 channels from 400 nm to 
850 nm for each station, there were about 180 equations available. 

By using the above method, R rs curves measured from waters of Monterey Bay 
and the Gulf of Mexico were inverted to determine the unknowns. The results and 
discussion are summarized in Chapter 7. 


*\ 
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CHAPTER 7 

RESULTS AND DISCUSSION OF THE INVERSE PROCESS 


For R n curves measured from waters of Monterey Bay and the Gulf of Mexico, 
the parameters a 0^(440), X and V, as well as the total absorption coefficients 
a(440), a(486) and a(550), were derived by minimizing the a.p.d. for each station. 
The variable estimates were obtained using Quattro Pro 5.00, by applying the 
quadratic feature. 

For R n measurements, the measured L sky may not be from the same part of the 
sky as that reflected off the sea-surface and entering the sensor when L u { 0 + ) is 
measured due to the sea-surface roughness. Cloudy conditions can worsen this 
mismatch. Therefore, instead of forcing r = 0.018, we let r and A in Eq. 39 be 
variables in the R rs inversion; thus r and A were also derived when the a.p.d. was 
minimized. For the 45 stations available, the overall average a.p.d. is 3.4%. 

From Eq. 27, a(440), a(486) and a(550) were also derived from the measured 
KJiav). As * 0.05 bja (Eqs. 23 and 29), 


K/av) 

l.0SD/0)(U20RJ ‘ 


(56) 


In the calculation of a by Eq. 56, two sets of D d { 0) were used. One was l/cos(/), and 
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the other was derived by forcing the ^(cv)-derived a(K) to approximate a w (\) + a p (X) 
+ a ? (X) for X > 600 nm. In this way, errors due to sea surface roughness (wave 
focusing of light) and ship shadow can be reduced, and DJO) can estimated for 
cloudy days when j is uncertain. 

7.1 Comparison of a(440), fl(486) and fl(550) Derived Using R n and K d Methods 

Figures 17a - 17c compare the results of the total absorption coefficients derived 
from the ^-inversion and those from K d methods (using derived at 440 nm, 

486 nm and 550 nm, respectively. From Figure 17a, it can be seen that K rj -derived 
o(440) is consistent with the ^-derived o(440), with r 2 = 0.94 (n=45), a slope equal 
to 1.03 and an average difference of 31 % for o(440) in the range of 6.03 m to 2.5 
m -i However, the average difference dropped to 19% for o(440) less than 0.5 m , 
perhaps because patchy data are more likely for turbid stations. When l/cos(/) was 
used to replace D d { 0) to derive o(X) from K d (\), the difference between the ^-derived 
and /^-derived o(440) was 73% for the whole o(440) range. This indicates that the 
method used to derive a(k) from K d (K) is preferable. These results also demonstrate 
that the method to obtain a(K) from R rs (K) inversion works very well even for this 
wide range of water types. 

The 31% difference between the results can be caused by the following factors: 
1) the errors in the measurements of L u , L sky , E d ( 0 + ) and E d (z), which will be 
transferred to R rs and K d , 2) simplifications in the model development;^) errors in fl. 
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a(440) from K d (m *) 


Figure 17a. Comparison of /?„-derived to A^-derix d a(440). 
For Details See Text. 


simulation; 4) water inconsistencies between the R n and K d measurements (temporal 
and spatial patchiness); and 5) method to obtain a from K d . With the consideration of 
these possible sources of error, a 31% difference seems rather small, and it can be 
claimed that we can not only qualitatively, but quite accurately derive the in-water 
absorption coefficient using remote-sensing techniques. 

Figures 17b and 17c compare <2(486) and <2(550) derived using the two methods. 
r 2 for <2(486) is 0.97 (slope 1.04) between the two sets of results, with an average 
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Figure 17b. Comparison of /{^-derived to /^-derived a(486). 

For Details See Text. 

error of 21%, while r 2 for a(550) is 0.97 (slope 0.87) with average error of 25%. The 
0.87 slope for a(550) means that the a(550) derived by R n inversion is consistently 
lower than that derived by K/av), with the implication that most of the 25% 
difference are systematic rather than random. One speculation is that simulated 
a^(550) might be low because of a large F value. This does not always happen, 
however; for many clear water situations, where phytoplankton particles are small and 
less "package effect" occurs, the simulated a*(550) is usually greater than the 
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0.1 1 

a(550) from K d (m *) 


Figure 17c. Comparison of i?„-derived to /^-derived a( 550). 

For Details See Text. 

measured o 0 (55O). Furthermore, for clear water, most of a(550) come from 0^(550) 

( *0.064 m' 1 ), so error in c^(550) simulations has only little influence on a(550), and 
it could not explain the 0.87 slope for clear water. Other possible sources of errors 
include the pure water absorption coefficient, the measurement of E d (z), and the 
possibility that the actual D/0) for 550 nm might be higher than the derived values 
since D d increases with scattering [Kirk 1991]; and there were relatively higher bja 
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ratio at 550 nm for many waters. However, how much higher D d (0, 550) should be is 
not clear. Further study is needed on this issue. 

7.2 Comparison of f?„-derived o(440), a(486) and <z(550) Using 
Simulated a* and Measured a p 

In order to see how the simulated a* curve performs compared to the measured 
a p curve, the R rs inversion process was also undertaken using the measured a p curve 
for each station. This time, it was assumed that we knew the curvature (shape) of a p 
for each station, but not the magnitude. The factor/was applied to the measured a p 
with no limitation set for its range. /was then derived by minimizing a.p.d.. Eq. 9 
was used for a g (\) because a d was contained in a p . The range for S g was set as 0.013 
< s < 0.017 nm' 1 , and the same range for Y as in section 6.3 was used. By the 
same process described in section 6.3,/ a/440), S g , X and Y were derived by 
minimizing a.p.d.. 

Figures 18a - 18c compare the results of ^-derived a(440), a(486) and 
a(550) values determined using the measured a p curves versus those determined using 
the simulated a „ curves. It can be seen that the results at all three wavelengths show 
close agreement. For the mathematical a * simulation (Eqs. 51 - 53), the r 2 values for 
a(440), a(486) and a(550) are 0.99 (n=48), while the average difference is 15.2% for 
fl(440), 10.3% for a(486) and 10.5% for a(550). For the empirical relationships 
(Eq. 55), the r 2 is 0.98 for a(440), and 0.99 for a(486) and a(550) (n=48), while the 
average difference is 23.2% for a(440), 18.6% for a(486) and 10.6% for a(550). 
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a(440) by a p curve (m *) 


Figure 18a. Comparison of /?„- derived a( 440) Using Simulated a* and 
Measured a p . 

The empirical a * simulation causes higher differences for the three absorption 
coefficients, which might be due to the simplification in the a * expression. From these 
results, we can say that the simulated a 4 curves work very well for the R n inversion 
process and can be used to replace the measured a p curve. 

In this process, the factors / for each station were also derived. For all 
stations, it varied from 0.18 to 2.34 with an average of 0.94, and more than 70% of 
the values fell in the range 0.7 to 1.3. Part of the / variation can be explained by 
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a(486) by a p curve (m 1 ) 


Figure 18b. Comparison of /?„-derived a(486) Using Simulated a ^ 
and Measured a p . 


patchiness and vertical structure, as the measured a p from one position may not 
represent the effective particle absorption coefficient of the upper water column, 
especially when a strong vertical structure exists. Additionally, /includes the errors 
introduced in the model development, variation of a w , and compensation among the 
parameters. 

For ideal situations, i.e. the water column is well mixed, correct values for a w , 
a g and scattering are used, and the ^-inversion process is well performed, / will be a 
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a(550) by a p curve (m *) 

Figure 18c. Comparison of /?„-derived a(550) Using Simulated a* 
and Measured a p . 

factor correcting the "beta factor" used in the a p calculation, since it is difficult to 
determine which "beta factor" is appropriate for our field samples (see Table 2 and 
Figure 12) as discussed in section 4.2. For offshore waters, pigment absorption 
coefficient is generally low. If the volume of filtered sea water is not enough, it is 
quite possible that the filer pad OD P is less than 0.10 (a p equivalent is -0.06 m 1 if 
the volume is 500 ml with GF/F 25 mm filter and using the "beta factor" of Bricaud 
and Stramski [1990]) where fi ^ vanes most. Thus, for better measurement of a p 
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derived ct^(440) (m *) 


Figure 19. Comparison of /J^-derived ^(440) to the Measured a*(440). 


regardless of the "beta factor" used, a large volume of water sample must be Filtered 
in order to reduce the "beta factor" uncertainties; this is consistent with a recent 
independent study by Patch et al. (in preparation). 


7.3 Comparison of ^-derived a* (440) to Measured a*(440) 


Using the simulations described in section 6.2, surface layer a 4 values 
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could be derived just from R n without prior knowledge of the "package effect" or 
the chlorophyll-specific absorption coefficients. Figure 19 compares the derived 
a*(440) to the measured 0^(440). The r 2 is 0.55 (n =36), but it becomes 0.92 (with 
46% difference) if/ is applied to the measured values, where the / values were 
derived as in section 7.2. As pointed out by Gordon and Clark [1980] and Gordon 
[1992], derived values from ^-inversion are optical averages of the upper water 
column, which can not be directly compared to in-water measurements made at a 
discrete position when the water is not homogeneous. Hence, /might be needed to 
compensate for possible differences due to water patchiness, vertical structure, and/or 
the "beta factor." These results indicate that at this stage, the <^(440) values derived 
by the ^-inversion method can be accurate to within 50% of the phytoplankton 
absorption at 440 nm in the upper water column. 
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CHAPTER 8 
APPLICATIONS 


8.1 Introduction 

Since an overflight can view a large area of the ocean in a short time period, 
applications of remote sensing via aircraft or satellite have become increasingly 
popular and important. These applications include surveys of pollutant flow, sediment 
transport, ocean circulation, estimation of [chi a] and [YS], modeling the water color 
at depth, and estimation of primary production, etc. In the following sections, some 
of these applications will be discussed, as well as the modeling of upwelling radiance 
of Tampa Bay measured from a low flying aircraft. 

8.2 Estimation of Chlorophyll Concentration 

The amount of phytoplankton in the ocean, often measured by the concentration 
of chlorophyll a ([chi a]), contributes significantly to regulation of the global climate 
system through its effect on the carbon cycle. It also indicates the trophic status of 
waters, and contributes to the conversion of light into heat. 

Since the 1970’s, efforts have been made to remotely estimate [chi a] in the 
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ocean. Empirical and semi-analytical algorithms have been reported for different 
waters and seasons [Gordon et al. 1980; Gordon and Morel 1983 and references cited 
there; Carder et al. 1991; Carder et al. 1994]. 

Applying the analytically derived <440) and the relationship between <440) 
and <675), [chi a] can be derived through 

[chi a] = -5 ? , ( 57 ) 

<( 675 ) 

where <(675) is the chlorophyll-specific absorption coefficient at 675 nm. 

In addition to the <440) value derived from /?„-inversion, the derivation of 
[chi a] depends on <(675) if we use Eq. 57. This value varies regionally and 
seasonally, but it has been shown to be more stable than <(440) as there is less 
influence due to the "package effect" at 675 nm [Carder et al. 1994], Thus, if we 
know the value of <(675), [chi a] can be derived using Eq. 57 with the R rs derived 
<440). 

It is necessary to keep in mind that /?„-derived [chi a] is an optically averaged 
value of the upper water column [Gordon and Clark 1980, Gordon 1992], Only when 
the water column is well mixed, can )?„-derived [chi a] be directly compared to 
measured [chi a ] from a discrete location. 

8.3 Estimation of the Gelbstoff Absorption Coefficient 


Generally, a d and a g have similar absorption spectra in the visible wavelengths, 
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so it is not easy to clearly separate a g from a dg , unless by active remote sensing such 
as laser induced gelbstoff fluorescence [Hoge et al. 1993], The scattering effects of 
detritus and gelbstoff are different, however, making it possible to estimate ^440) 
from the X or R n values. Then a g ( 440) can be separated from the ^-derived a dg ( 440) 
as discussed in sections 6.3 and 7.1. 

V 

Figure 20 shows the relationship between measured a d { 440) and the derived X 
values. In the log-log format, r 2 between a d { 440) and X is 0.89 (n = 39); and the best 
regression fit is 

a 1 440) = 61.44X 1 31 . (58) 

Thus, with the 0^(440) and X values derived by R rs in section 7.1, 

^(440) = 0^(440) - 61.44X* 31 . 

Eq. 58 might be strongly driven by the sediments from the Mississippi River. 
For other regions, different regression results may occur. Also, there may be a 
difference between the "beta factor" of detritus and that of phytoplankton [Bricaud 
and Stramski 1990, Nelson and Robertson 1993]. More data sets and better 

measurements are expected to improve Eq. 58. 

As ^(440) = a g \ 440) * [YS], it is possible to estimate the gelbstoff amount if 
we know the gelbstoff-specific absorption coefficient, a g (440). Unfortunately, as with 
the chlorophyll-specific absorption coefficient, a g (440) varies for different sources of 
gelbstoff [Carder et al. 1989, Hawes 1992], making it difficult to estimate [YS]. How 
fl ^*(440) varies is beyond the scope of this research, although we can choose an 
averaged £^*(440) value for the calculation of [YS]. 
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Figure 20. ^-derived X vs. Measured a/ 440). 


8.4 Algorithm for the Absorption Coefficient at 490 nm, a(490) 

Optical water types can be classified by the attenuation or absorption coefficient 
at 490 nm [Jerlov 1976; Austin and Petzold 1981], and 1/^490) yields a measure of 
the penetration depth of solar light [Gordon and McCluney 1975] in meters. By 
analyzing the diffuse attenuation coefficient #(490) (which is proportional to #^490)) 
and the in- water upwelling radiance ratio at 443 and 550 nm, Austin and Petzold 
[1981] developed the algorithm 
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£(490) = 0.0883 


L„(443) 
v 1.(550) , 


k-1.491 


► 0.022 , 


(59) 


where the 0.022 is the attenuation coefficient of pure sea water at 490 nm [Austin and 
Petzold 1981] in units of m' 1 . 

Applying this algorithm to the West Florida Shelf data and other stations, 

Carder et al. [1992] found that the correlation between <3(490) and the ratio 
/? rj (442)//? ra (550) is less than that between <2(490) and /?„(520)/£„(560) (Figure 21a). 
This may partly be explained by the influence of bottom reflectance and other 
components that do not co-vary with chlorophyll for "case 2" waters. 

For regions that include shallow, coastal and "case 2" waters, an algorithm 
similar to Eq. 59 is presented for the calculation of <2(490) based on the study of 45 
stations. These stations include blue-sky and cloudy sky conditions. Alternatively, 
a(490) could be analytically derived by using high-resolution R rs values as discussed 
in sections 6.3 and 7.1. 

For the algorithm development, /^-derived <2(490) (= o(486), using derived 
D/0)) is considered as the true value. Figure 21a shows the relation between the 
reflectance ratio and /^-derived <2(490). For <2(490) ranged from 0.03 to 1.5 nv\ it is 
found that r 2 is 0.96 (n=45) between ln(a(490)) and ln(7?„(520)//? ra (560)), while r 2 is 
0.89 between ln(c(490)) and ln(i?„(442)/2?„(550)). Thus, for the best fit, an algorithm 


for <2(490) is 



ratio derived a(490) ( 
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a(490) = 0.19 


RJ 520 )\ 
**( 560 ), 


- 3.11 


( 60 ) 


with average error of 32.5%. However, if the ratio /f„(442)/i? rj (550) is used instead, 
the best fit generates 


<j(490) = 0.15 


(R„{ 442) \ 
,**( 550 ), 


- 1.37 


( 61 ) 


with average error of 49.6%. Remember that the error was only 21% when the 
hyperspectral R„-inversion method (section 7.1) was employed. Figure 21b shows the 
results from Eq. 60 and Eq. 61. 


8.5 Modeling the Water Color at Depth 

Combining Eqs. 27 and 44, one can derive 

B& - El o>- 108 ^ . (62) 

Since DJd 0) * 1/cos (/) [Gordon 1989a], E d (z) can be modeled if E d ( 0 ) and a are 
known. From sections 6.3 and 7.1, a can be derived from R rs inversion, and E d (0) 
can be calculated by models [Gregg and Carder 1990, Bishop and Rossow 1991]. As 
examples, Figures 22a - 22d show the modeled versus measured E d (z) for station 
GO10 at wavelengths 440, 486, 520 and 550 nm, respectively, with £^(0") 
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Figure 22d. Measured vs. Modeled ^(550) of Station GOIO 
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calculated by the model of Gregg and Carder [1990]. It can be seen that the model 
gives very close subsurface irradiance (about 10% difference) compared to the 
measured values; and the calculated irradiance profiles were consistent with the 
measured for the upper water column. 

8.6 Estimation of Primary Production 

In the past decades, primary production models have been based on pigment 
concentration, either by light available (P-I relationship of Platt et al. [1991]) or light 
absorbed (P-AQ relationship of Bidigare et al. [1992]), using spectral or non-spectral 
expressions. Recent popular models use the light absorbed approach, i.e.: chlorophyll 
a concentration multiplied by the irradiance and two factors: averaged chlorophyll- 
specific absorption coefficient, and the quantum yield for carbon fixation. Thus the 
photosynthesis rate at depth z is [Kishino et al. 1986, Smith et al. 1989, Cullen 1990, 
Morel 1991, Marra et al. 1992, Zaneveld et al. 1993]: 

PP(z) = 4 *z)[chl a]alPAR(z) , (63) 

where 4>(z) is the quantum yield for carbon fixation, and [chi a] is the chlorophyll a 
concentration (mg/m 3 ). is the spectrally-averaged chlorophyll-specific absorption 
coefficient over 400 - 700 nm (m 2 (mg chi a)' 1 ): 
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_ J 400 [chi a\ 

°+ r°°E 0 \z,X)dX 

J 400 

and PAR{z) is the photosynthetically available irradiance, computed by 


(64) 


PAR(z) = r°°E 0 \z,k)dX , 

J 400 


(65) 


where E 0 q (z,\) in Equations 64 and 65 is the spectral scalar irradiance in 
quanta/ m 2 /nm/s. From Sathyendranath and Platt [1988], 


E 0 \z,X) * E o \0-,X)e-*^ . (66) 

Quantum yield 4>, a measure of the phytoplankton growth rate [Platt 1986], 
varies with light intensity, physiological status, temperature and nutrient stress of the 
phytoplankton population [Langdon 1988, Cleveland et al. 1989, Smith et al. 1989], 
We choose the empirical formula suggested by Kiefer and Mitchell [1983] to express 

how <f> changes with PAR, 


Jhl \ - A + (67] 

PAR^PAR(z) ’ 

where 4> m is the maximum quantum yield, PAR* ( K + in Smith et al. 1989) is the PAR 
value at which <j> = <t>J 2. 

Therefore, combined with the photoinhibition expression suggested by Platt et 
al. [1980], primary production at depth z can be expressed as 
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with 


PP(Z) - AP(z) *- PARt f' Fm « , 
PARq+PAR(z) 


(68) 


AP(z) = y [C]a;PAR(z) , (69 

where v in Eq. 68 is a parameter to describe photoinhibition. 

Eq. 69 accounts for the absorbed photons by the phytoplankton pigments, in 
which 7 determines the fraction of chlorophyll a relative to the total pigment and 
phaeo-pigment concentration [CJ. 

Traditional estimates of primary production based upon remote measurements 
are accomplished using Eqs. 68 and 69, i.e. [ C] is estimated first from remotely 
measured signals [Vargo et al. 1987, Platt et al. 1991, Balch et al. 1992] using, for 
example, the CZCS algorithm: 


IQ 


= A 


Z- w (443)>^ 


L( 550) 


* 0.95^4, 


'*„( 443) 
*„(550) 


y*2 


( 70 ) 


where A, = 1.13 and A 2 = -1.71 [Gordon et al. 1983]. The 0.95 value comes from 
Ej(443)/E d (550) = 0.95. With this estimated [Q, the diffuse attenuation coefficient of 
irradiance can be approximated by the empirical relationship suggested by Morel 
[1988]: 


K/k) = K w (\)+d(k)[C]‘ w , 


( 71 ) 


in which values for K w , d and e are found in tables calculated by Morel [1988]. 

With the measured or calculated E 0 q { O') and K d , PAR at any depth can then be 
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estimated. However 7 , v, <f> m , PAR + and a* must be estimated for the calculation of 
PP. 7 varies from 0.3 to 0.9, with an average of 0.75 [Morel and Berton 1989, Balch 
et al. 1992], the value used for our calculation. 4> m , v and PAR * vary with 
phytoplankton physiological status, which at this point cannot be derived based upon 
remotely sensed measurements and must be estimated from other information. 4> m has 
been reported to range from 0.03 mol C (Ein absorbed)' to 0.1 mol C 
(Ein absorbed ) -1 [Bannister and Weidemann 1984, Smith et al. 1989, Morel 1991]. As 
a kind of average for productive waters, 4> m is assumed to equal 0.074 mol C 
(Ein absorbed) 1 , a value suggested by Cullen [1990]. There is little literature 
information, however, about the values of v and PAR *. Platt et al. [1980] found v 
varied from 0 to 0.0028 (W/m 2 )" 1 . If we choose 0.0028 (W/m 2 )’ 1 for v, its equivalent 
value is v * 0.01 (Ein/mVday)' 1 . Kiefer and Mitchell [1983] found that PAR* equals 
about 10 Ein/m 2 /day, a value assumed for the waters in this study. 

The only unknown remaining is a‘. This value varies with phytoplankton 
condition (a* (A)) and the light environment [Morel 1978, Kishino et al. 1986]. If 
only remote-sensing information is available, a ^ must be estimated from other data, 
or an empirical value has to be picked for a/. In many studies, an a* value of 
0.015 (mg chi a/m 2 )- 1 has been used [Bannister 1974, Dubinsky et al. 1984, Smith et 
al. 1989 and Marra et al. 1992] and assumed to be vertically constant. 

Using the CZCS algorithm, [C] is found to be accurate to about a factor of 2 
[Gordon and Morel 1983], and a/(440) can vary by a factor of 4 [Morel and Bricaud 
1981, Bricaud et al. 1988, Laws et al. 1990, Stramski and Morel 1990, Carder et al. 
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1991] due to pigment composition, "package effects," and the color of the light field. 
That means the combined variation of [C] and a / for AP could be off on a global 
scale by a factor of 8 when other terms are known, although the range of uncertainty 
is probably significantly less than 8 here as we are not in subtropical waters (e.g. see 
Laws et al. 1990). 

Eq. 63 is actually a simplified version of a complete spectral expression as 
noted by Sathyendranath et al. [1989b], Morel [1991] and Platt et al. [1991]; i.e. 

PP(z) = 4>(z)f™a^\)E 0 \z,k)d\ , ( 72 ) 

and the absorbed photons are 

AP(z) = f 10C a (k)E o \z,\)dk . ( 73 ) 

J 400 v 

However, those full spectral models were based on the pigment or chlorophyll a 
concentration, similar to Eq. 63. If a 0 (X) and o(X) can be obtained directly and 
analytically from remotely measured signals (e.g. remote-sensing reflectance), instead 
of using the pigment concentration based model and the empirical relationships to 
obtain [C] and the attenuation coefficient for AP(z), the problems involved in the 
estimation of [C\ and choosing values for 7 and a / will be avoided, and the accuracy 
of estimating AP and PP would be improved when O'), <i> m , v and PAR + are 
certain. 

From the discussions in Chapters 6 and 7, we know a^(X) and a(X) can be 
derived solely from measured ft ra (X). Thus, scalar irradiance at depth z can be 
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calculated through Eq. 66 with K d ~ 1.08 D d (0) a [Gordon 1989a]. E o 9 (0) can be 
determined from models of Gregg and Carder [1990] or Bishop and Rossow [1991]. 
Therefore, when 4> m , PAR+ and v are known, PP for any depth can be calculated 
through Eqs. 68 and 73. 

For data collected from the Marine Light-Mixed Layer (ML-ML) study 
(21°W/59°N) of May 1991 ( PAR and primary production measurements are in Marra 
et al. [1994]), two sets of calculations were made. One set is derived using Eqs. 68 
and 69 (referred to as the "pigment" method hereafter), as discussed at the beginning 
of this section. The other set is derived by Eqs. 68 and 73 (referred to as the 
"absorption" method hereafter), where the absorption coefficients of the pigments and 
the total absorption coefficients are analytically derived from measured remote-sensing 

reflectance. 

In both calculations, measured PAR(0) for each station was used for the 
calculation of PAR(z) and comparisons of the two methods, because most of the 
stations were taken during cloudy weather. D d ( 0) is approximated as 1.2 [Platt et al. 
1991] in the "absorption" method as most of the stations were under clouds. The 
calculated results for each day at each depth were compared with the measured 
primary production values, and the water-column-integral production. 

Comparisons of calculated to measured PAR and calculated to measured PP are 
presented in Figures 23a - 23d and Figures 24a - 24d, respectively. All data are 
summarized in Table 7 and all PP data in Figure 25. 

It can be seen that the "absorption" method yields close estimates of PAR(z), 
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while PAR(z) values calculated by the "pigment" method were consistently higher, 
particularly at depth. One reason for this is that the CZCS pigment algorithm 
estimated [chi a ] as much as a factor of 5 lower than the measured surface values for 
these waters. This had the consequence that K d values calculated from these [C] 

([chi a]/0.75) were small. This difference could be the result of an incorrect CZCS 
algorithm for that environment [e.g. Balch et al. 1989, Mitchell and Holm-Hansen 
1991], or less likely there were substantial errors or discrepancies in the 
measurements of R n or "sea truth" [chi a]. Additionally, the empirical relationship 
between K d and [C] (Eq. 71) might not hold for these waters. The same field data, 
however, are used for both methods; so models and algorithms will be largely 
responsible for differences between the methods. 

Primary production values calculated using the "absorption" method were highly 
comparable to the measured values, whereas the "pigment" method significantly 
underestimated PP, particularly in the surface waters. The average difference of the 
water-column-integrated PP is 20% between the measured and the "absorption" 
method, while it is 61% between the measured and the "pigment" method (Table 7). 
The r 2 is 0.95 (n=24) for a linear regression between the values calculated by the 
"absorption" method and that of the measurements, with a slope of 1.26 and + 32% 
difference; whereas the r 2 is 0.85 (n=24) between the values calculated by the 
"pigment" method and that of the measurements, with a slope of 0.34 and - 78% 
difference (see Figure 25). These results indicate that there is a factor of 3 
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Measured vs. Modeled PAR of May 17 in ML-ML. 



Figure 23b. Measured vs. Modeled PAR of May 20 in ML-ML. 
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. Measured vs. Modeled PAR of May 22 in ML-ML. 



Figure 23d. Measured vs. Modeled PAR of May 24 in ML-ML. 
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Table 7. Results about PP Calculation of ML-ML, May 1991. 


Day May 17 

May 20 

May 22 

May 24 

PAR(0) (Ein/m 2 /day) 

38.27 

16.25 

65.72 

28.73 

R„( 443)/R„(550) 

1.2 

1.6 

1.8 

1.8 

surface [chi a] (mg/m 3 ) 

2.9 

1.3 

1.5 

1.0 

CZCS [chi a] (mg/m 3 ) 

0.67 

0.43 

0.36 

0.34 

Mea. integral PP 
(m mol C/m 2 /day) 

190 

98 

89 

99 

Est. integral PP: "absorption" 
(m mol C/m 2 /day) 

235 

94 

105 

89 

Est. integral PP: "pigment" 
(m mol C/m 2 /day) 

77 

34 

62 

40 


improvement in the accuracy of PP calculation from remote sensing for ML-ML 

waters using the "absorption" method. 

For May 17, the calculated PP values by the "absorption" method in the 
euphotic zone were generally greater than the measured values. Factors that can 
account for these differences include possible errors in measurement and R rs 
inversion, grazing effects on the measured values, and either the photoinhibition 
factor used was smaller than the "real" situation, or 4> m and PAR * used were higher 
than the "actual" values. 

For May 20 and May 24, the PP values calculated by the "absorption" 
method were very close to those measured, which suggests that our assumed values of 
<t> m , PAR * and v for these two days were close to the realistic situation. 

The greatest differences between measured and calculated values from both the 
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Figure 24a. Measured vs. Modeled PP of May 17 in ML-ML. 




Figure 24b. Measured vs. Modeled PP of May 20 in ML-ML. 
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Figure 24c. Measured vs. Modeled PP of May 22 in ML-ML. 
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Figure 24d. Measured vs. Modeled PP of May 24 in ML-ML. 
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Figure 25. Measured vs. Modeled PP of ML-ML, May 1991. 


"pigment" and "absorption" methods occurred on the only sunny day of the cruise, 
May 22, where measured surface layer production was smaller than in the second 
layer. Other than measurement errors, this is an indication of increased 
photoinhibition, as suggested by the calculated values of the "absorption" method. 

The calculated PP values by the "absorption" method, however, are greater than the 
measured values for the first four depths (~40%). Possible reasons for this difference 
might be an overestimation of <f> m or an underestimation of v for the phytoplankton in 
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this environment, since <t>„ and v might be a function of light history. The other 3 
cruise days were overcast with an expected adaptation of the plants to a low light 
environment. The sudden exposure to the bright light probably caused extreme 
photoinhibition in the dark-adapted phytoplankton population with the resultant 
overestimation of PP by the "absorption method. 

The PP values calculated by the "pigment" method for the sunny day, however, 
do not show a photoinhibition response in the surface waters (Figure 24c). One reason 
for this is that ~aj not only varies with the chlorophyll-specific absorption coefficient, 
but it also varies with the light field (Eq. 64). An a/ value for surface water is not 
necessarily appropriate for the whole water column, even for well-mixed oceans 
[Kishino et al. 1986]. So, the ~a* value is one of the major sources of error in PP 
calculations by the traditional method. 

For the waters studied, PP values calculated by the "pigment" method were 
significantly lower than the measured rates. Since PP values calculated by the 
"absorption" method were close to measured values, our value of <j> m , 0.074 mol C 
(Ein absorbed)' 1 , is apparently close to the actual value. Therefore, the differences 
between the "pigment" method calculated and measured PP are largely due to the 
estimation of [C] and the product of 7 {C]a/. However, when only remote-sensing 
data are available, it is difficult to know which are the "correct" a<, and 7 [C] values, 
as they depend upon empirical parameters. What makes things interesting is that a 
factor of 2 increase in both [C] and 7/ will make calculated PAR and PP both close 
to the measured values. Table 7 pigment data, however, suggest that the CZCS 
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algorithm accounts for most of the error. This is consistent with the speculation of 
Platt et al. [1988] that determination of biomass by remote sensing dominates the 
error in primary production estimation. 

Comparing Eqs. 69 and 73, the main difference between the "absorption" and 
"pigment" methods is how AP(z ), the absorbed energy by phytoplankton pigments, is 
obtained from remote-sensing data. In the traditional method, calculation of AP 
depends heavily on 4 parameters: 7 , A,, A 2 , and a^. When only remote-sensing 
signals exist, it is difficult to verify which values should be used for each of these 
parameters, although each of them could be tuned empirically to specific regions and 
specific seasons [Platt et al. 1991]. Since there are errors associated with each 
number, the cumulative error in AP could then be very high even if PAR is correct. 

In the AP calculation by the "absorption" method, however, the phytoplankton 
absorption and the total absorption are analytically derived from measured R n . Thus, 
most of the error comes from the /^-inversion process. It is believed that absorption 
coefficients derived from R rs inversion have an accuracy better than 50% for "case 1" 
and "case 2" waters (Chapter 7). This means that the accuracy in the AP calculation 
by the "absorption" method could be within 50% compared to the errors as great as 
400% derived using the traditional method. In this "absorption" method suggested 
here, the mid-step, i.e. deriving [chi a] and estimating a value for a*, is eliminated. 
The result is highly improved accuracy in deriving AP and PP. The remaining 
challenge for the calculation of PP through remote sensing is how to remotely obtain 
accurate estimation of the physiological parameters [Balch et al. 1992] such as 4> m , 
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PAR^ and v, that could perhaps be based on light history, chlorophyll fluorescence 
efficiency [Chamberlin et al. 1990], and surface temperature anomalies [Kamykowski 

and Zentara 1986]. 

In conclusion of this section: 

a) Based on the results presented here, the a 9 simulation works very well in 
modeling PAR(z) from PAR( 0) and measured R ra . Combined with parameters 
regarding photosynthesis, the PP values calculated by the "absorption" method were 
close to the measured ones with r 2 = 0.95 and slope = 1.26; which is a factor of 3 
improvement in PP estimation accuracy over the traditional method. These results 
also indicate that the combined PP model, 


PP(z) = AP(z) 


PAR.+PAR(z) 


-vPAR(i) 


(74) 


works well for the ML-ML waters. 

b) It is not necessary to know [chi a\ and a * for the calculation of PP. In the 
traditional method, the estimation of [chi a] (or [Q) and a,* separately is an inherent 
disadvantage in calculating PP based upon satellite or aircraft remote sensing of ocean 
color. From the remote-sensing point of view with regard to remote estimation of PP, 
what is really needed is the absorbed energy by phytoplankton and the physiological 
parameters regarding the photosynthesis by phytoplankton. So, it is preferable to shift 
the primary production model from the pigment-concentration based to the pigment- 


absorption based. 
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8.7 Modeling Total Radiance (L,) Measured from a Low-flying Aircraft (Blimp 
Shamu) 

In January 1993, the total radiance, L„ of Tampa Bay was measured from the 
SeaWorld blimp "Shamu". These L, measurements covered waters over sandy and 
grassy bottoms and dark, gelbstoff-rich waters. Although values of E d were measured 
before and after the flight, no E d or proper L sky data could be collected when L t was 
measured, so no precise R n could be derived using standard methods. However, L sky 
from part of the sky was measured, and its spectral shape might be considered 
approximately correct. Then, with the method discussed in Chapter 7, the parameters 
r and A were used to estimate the sky radiance which enters the sensor. Thus the 
measured L , could be modeled, and the in-water absorption coefficient and the water 
depth for optically shallow water could be derived. 

The total radiance entering a sensor in clear air is [Gordon and Wang 1994] 

I, = L u ( 0>' ( ^ To) + L r + L a , < 75) 

where L R is the contribution due to Rayleigh scattering in the atmosphere, L A is the 
contribution due to aerosol scattering, and exp {-t a - r c ) determines the attenuation of 
atmosphere to the upwelling radiance. 

Since the altitude of the blimp was only about 50 m over the water surface, 
exp(-T< - t g ) is * 1 and the small contributions due to L R and L A are combined into 
rL sky and A E d of Eq. 40, so 
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Figure 26a. Measured vs. Modeled L, of Station SH04. 



Figure 26b. Measured vs. Modeled L, of Station SH34. 
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Figure 26c. Measured vs. Modeled L, of Station SH35. 


L. * (R „* A ')~!r* rL *y ■ 

k g 


( 76 ) 


L g was estimated by interpolation between the pre- and post-flight 
measurements of L a . With the a* simulation discussed in section 6.2, L, of waters 
over sandy bottoms (SH04), grassy bottoms (SH35), and dark waters (SH34) were 
modeled. Since R n and L, are interchangeable when other terms are certain, R n curves 
of above situations were also modeled at the same time. 

Figures 26a - 26c show the measured and modeled L„ while Figures 
27a - 27c show the measured and modeled R n . It can be seen that excellent model 


results were achieved for these waters. The bottom albedo of the sandy bottoms came 
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Figure 27b 


Measured vs. Modeled R rs of Station SH34. 
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Figure 27c. Measured vs. Modeled R„ of Station SH35. 


from Figure 10, while the bottom albedo of the grassy bottoms came from reflectance 
of vegetation [Tucker and Miller 1977]. 


bottom albedo 
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CHAPTER 9 

CONCLUSIONS OF PART B 


9.1 The phytoplankton absorption coefficient, o*(X), can be simulated using 
the following expressions: 

1) mathematical functions: 

400 i ^ 570, a + (A) = a^e 



570 < X < 656, a^\) 


«*( 570) 


,,(656) - ^5 70) 5?0)j 

656 - 570 


(78) 


and 


(x - x 2 ) 2 

. 2a\ 

656 < X < 700, a^X) = a^e 


(79) 


2) empirical relationships: 

a*(X) = [fl 0 (X) + a 1 (X)ln(fl^(440))]^(440) . (80) 

The mathematical expression involves 6 parameters, but only 4 of them (a*,, F, 
a +2 and a 2 ) vary significantly among different waters. 3 of the 4 parameters, F, o 2 and 
a^/a^i could be estimated from the value of a ^ through 
= 0.86 + 0.161n(o^), 
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F = 2.89exp[-0.505tanh[0.561n(^/0.043)]], 

a 2 = 14.17 + 0.91n(aJ. 

a 0 (k) and a,(X) of Eq. 80 for each wavelength were empirically derived 
(Table 6). So, using either method, a*(\) curve can be simulated when a*(440) is 
known. Comparing the two methods, mathematical functions create a smooth a* 
curve, but with more potential to adjust the curve shape. The empirical relationships, 
however, provide empirical spectral finesse to the curves. 

9.2 With the above ^(X) expressions, the inverse problem of R rs can be 
reasonably well solved. For waters of Monterey Bay and the Gulf of Mexico, the 
/?„-derived total absorption coefficients at 440 nm, 486 nm and 550 nm were 
consistent with the values derived from K d . When using the mathematical simulation, 
the average difference is 31% for o(440), 21% for a(486) and 25% for <a(550) for 
fl(440) ranged from 0.03 nr 1 to 2.5 m' 1 . These results were also consistent with the 
values derived using measured a p curves. These results demonstrate that the fl^(X) 
simulations can be used in remote-sensing applications without prior knowledge of the 
in-water optical properties, such as chlorophyll-specific absorption coefficient. 

9.3 The ranges for the 5 variables (a^,, «*(440), S*. X, and Y) are close to 
realistic situations. 

9.4 With the ^(X) simulation, the estimation of primary production in the 
euphotic zone can be carried out just from measured R rs and surface PAR. The 
problems associated with the traditional method, i.e. estimating [chi a] and using a 
chlorophyll-specific absorption coefficient, are avoided. For data collected in the ML- 
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ML May 1991 waters, r 2 is 0.95 between the calculated and measured PP values, 
with a slope of 1.26 and 32% average error. This is a factor of 3 improvement in the 
accuracy of PP estimation by using the suggested method over the traditional method 
for those waters. However, challenges remain to remotely estimate the physiological 
parameters since they might be functions of nutrient stress, life stage, light history, 

etc. 
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CHAPTER 10 

SUMMARY AND FUTURE WORK 


10.1 Summary 

a. A remote-sensing reflectance (/?„) model is developed, including terms for 
the contributions of bottom reflectance, gelbstoff fluorescence, and water Raman 
scattering. 

b. The remote-sensing reflectance model is tested for a wide range of water 
types, and an average error of 3.4% is obtained between the measured and modeled 
remote-sensing reflectance curves. This result indicates that the forward problem, to 
interpret the water color or remote-sensing reflectance in terms of in-water 
constituents, is well solved. 

c. A six parameter model is developed to simulate the spectral absorption 
coefficient of phytoplankton pigments (fl 0 (X)) with an average error of about 11%. It 
can be used in the remote-sensing reflectance modeling and inversion. 

d. Using the ^(X) model and the ^-inversion methodology developed above, 
in-water absorption and scattering coefficients can be analytically derived solely from 
the remotely measured signals. For o(440) ranging from 0.03 to 2.5 m' 1 , the remote- 
sensing-derived a{ 440) values are within 31% of the in-water measured values. This 
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indicates that the inverse problem, to analytically derive in-water optical components 
from remotely measured signals alone, is reasonably well solved. 

e. A new approach is suggested for the estimation of primary production based 
on remote-sensing measurements. In this new approach, the primary production model 
is based on the phytoplankton pigment absorption coefficient, instead of the 
phytoplankton pigment concentration. 

f. Using the new approach for a study in the North Atlantic during May 1991, 
the estimated primary production at depth in the euphotic zone based on remote 
sensing was within 32% of the measured values. This result is about a factor of 3 
improvement in the estimation accuracy over a traditional method, which is based on 
the pigment concentration. 

10.2 Future Work 

a. Validation for values of X and Y is not complete, which requires independent 
in-water measurements of the values and distributions of the spectral volume 
scattering function and radiance field. Hopefully, this will be solved in the coming 
years with the development of new, sophisticated instrumentation. 

b. Further coupled in-water and remote-sensing measurements for a wide range 
of waters are necessary in order to improve and expand the usefulness of the model 
and simulations suggested in this study. 
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c. Part of the errors in the component derivation come from the internal 
compensations among the parameters. In order to improve the accuracy of derived 
in-water components from remotely sensed signals alone, theoretical and experimental 
studies about the parameters need be carried out. For example, how great is the 
compensation among the parameters? What should be the right ranges for Y and V 
How to determine their ranges? etc. 

d. Methods to accurately derive a(k) from K d (\) need further study. 

e. Coupled investigations between primary production and remote-sensing 
reflectance measurements need to be carried out widely to generate more data to test 
and improve the suggested new approach. 

f. For the remote estimation of primary production, methods must be pursued to 
remotely or empirically estimate changes of the physiological parameters, perhaps by 
their co-variance with some remotely measured variables such as sea-surface- 
temperature anomalies, wind-stress history, and light history, etc. 
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APPENDIX 1. EXACT SOLUTION OF NADIR R„(a) BY THE RADIATIVE 
TRANSFER EQUATION 


The radiative transfer equation for nadir radiance L u {z) is 


dLtf 

dz 


= ~cL(z) + L' u 


( 1 ) 


where L ’ is 


L* = f f}(a)L(dWj)du' . < 2 

If we separate the radiance field L into two parts: L d for the radiance in the 
downwelling field, and, L u for the radiance in the upwelling field, then Eq. 2 can be 
re-written as 


Define the light-averaged-backward-VSF and VSF-averaged-upwelling- 
radiance L u similarly to Zaneveld [1982, 1994], respectively: 

f pCaJL/Qfya' 

p _ J2it 

L i (0 / )cos(0V w/ 

f P(a)L u (Q)d«>' 

Y _ » -2 7i 

“ ~ C / 9 

f_ u 


then Eq. 1 becomes 
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dL(z) — _ 

—j— = ~cL u (z) + p b E d + L u (z)b f . 


(5) 


If we define the diffuse-attenuation-coefficient for nadir-viewed radiance L as 


k L (z) = - 


dL(z ) 


L u (z)dz ’ 

then from the R n definition, the nadir in-water remote-sensing reflectance /? rj (0 ) is 


( 6 ) 


**«n s 


c+k L (0~)-ebj 


(7) 


where 


€ 



( 8 ) 
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APPENDIX 2. DIFFUSE ATTENUATION COEFFICIENT OF AN EXTENDED 
LAMBERTIAN RADIANCE SOURCE 




Figure 28. Schematic light field of an extended Lambertian source. 


Figure 28 illustrates the light field illuminated by an extended Lambertian 
source. L g is the source radiance and it is the same for all directions. Define & as the 
diffuse attenuation coefficient of radiance, then the radiance at distance z in direction 
0 is 


fa 

L(0»z) - l/"® 


( 1 ) 


and the irradiance at that distance is 
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_K fe 

EJ,z) = 2nL g j*e cos(0) cos(0)sin(0)J0 . 
Notice that the source irradiance is 


EJP) = it L 0 , 


then. 


E/P) 


JL fa 

= 2( l e COi(0) cos(0)sin(0)d0 . 

JO 


From Gordon [1989a], we have 


( 2 ) 


(3) 


(4) 


E / z ) _ - 1 . 08 Dj( 0 )icz (5) 

^0) " 

and for a Lambertian distribution, £> d (0) = 2, so 

Jt fa 

e - 216 « * 2 f 2 e cos(6) cos(0)sin(0)<f0 . (6) 

Jo 

Define k = £k, kz = t, then for r in the range of 0.5 to 2.0, f varies from 1.37 
to 1.62 for Eq. 6 to hold. So, in a general case, f = 1.5 is used as an average. 

Also, from Eq. 36, for bottom albedo p = 0.5, r = 0.5 means R rs b * 0.022 sr' 
\ and t = 2 means R a b * 0.0004 sr' 1 . If the average remote-sensing reflectance from 
the water column is 0.005 sr' 1 , r = 0.5 means R rs b may contribute -99% of the total 
signal, and r = 2 means R n b contribute less than 10% of the total signal, which 
indicates 0.5 to 2 of r is the general range where remote-sensing reflectance from the 


bottom could be verified. 
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APPENDIX 3. REMOTE-SENSING REFLECTANCE FROM GELBSTOFF 
FLUORESCENCE (Rj) AND 
WATER RAMAN SCATTERING (/?„*) 


For z positive downward from the surface (Figure 5), with the consideration of 
isotropic (i u , to first order the inelastic radiance (L u ie ) in the direction 0 and the 
upwelling irradiance (E uje ) at depth z due to the depth interval dz are simplified to 

dL,Jz$M - «) 

and 


dE uU (z,X) = -2n f*dL u Jz,6, X)cos(Q)sin(e)d6 
= [ flJX x ,X)E o (z,X x )dX x dz, 

J X x 

where EJz.XJ is the scalar irradiance at depth z, and E 0 (z,\) = E^z.XJ + 

« D d (l+2R(XJ)E <i (z > XJ. Considering that R(k x ) is small (<0.05) and DJD d * 2 
[Gordon et al. 1975] is independent of depth, the subsurface irradiance due to the 
inelastic scattering for a deep water column is 


E„ . (0-A) . 2 it / 


P te (M)Ej(0-,XJ jv 

2k(X)+k(X x ) 


( 3 ) 


Defining Q ie as the "Q factor" for the inelastic scattering field, the subsurface 
upwelling radiance due to inelastic scattering is 
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2k(X)+k(X x ) 


dX. 


The inelastic total scattering coefficient (m'/nm) is defined as 


(4) 


»K**>*) = f 4n P u (<*>X x ,X)do). (5) 

Since 0 u (a,\ x ,\) is considered isotropic, then 

i|r(X x ,X) = 4nP ie (X x ,X). (6) 

According to the definition of remote-sensing reflectance, with Eq. 4 and Eq. 6, we 
have 


u« 


r 2 r 

2Q u n 2 J 1 *[2k(X)+k(X^ 0-,X) x ' 


(7) 


For gelbstoff fluorescence, defining rjfXJ as the quantum efficiency for the emission 
band excited by \ x , then [Gordon 1989, Carder and Steward 1985] 


n(.K) - L 


X jKjvjj 

* K 


dX. 


( 8 ) 


^(X,,X) can be characterized by a log-normal curve [Hawes et al. 1992], so 




XA 



(9) 


in which 



where \ f , s and a may vary with the type of gelbstoff and 

In general, b b < < a for most oceanic waters [Morel and Prieur 1977], so k is 
close to a. And, based on the calculation for chlorophyll a fluorescence made by 
Gordon [1989], the Q ie factor for inelastic scattering is ~3.7. Then combining Eq. 7 
and Eq. 9 with ( * 0.98, n w ~ 1.34, the remote-sensing reflectance due to gelbstoff 

fluorescence can be reduced to 


- o.onf k MK) 






X [2a(\)+aa )]£/ 0 -,X) A 


d\ x . 


(ID 


Unlike broad-band (- 100 nm) gelbstoff fluorescence, the water Raman 


emission has a half-band width of about 20 nm [Collins et al. 1984]. Omitting this 
band width, i.e. assuming a narrow Raman emission, the inelastic scattering 
coefficient for water Raman can be related to Raman scattering coefficient as 

\p(K x ,X)d\ = b*(K) 

and from Eq. 7, with k * a, the remote-sensing reflectance for water Raman is 


b R (k x )EJ0-,\ x ) 


K*(X) « 0.072 


[2a(X)+fl(X I )]£ <i (0-,X) 


( 12 ) 




